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INTRODUCTION 




Several current and anticipated trenda threaten to limit the uaefulneaa 
of future NASA remote aenaing miaaions if we continue to employ current 
data handling methods. Forecasts of instrumentation trends show an in- 
crease in sensor resolution from 80 meters for Landsat 1 to IS meters 
for the Operational Earth Resources Satellite (OERS) planned for launch 
in 1990 (Ref 1).* In c^ddition to the increased resolution, the OERS 
design currently calls for 20 different bands. This translates to an 
increase in total data lates for earth sensing missions from 1 million 
bits per second (Mbps) to over 100 Mbps, an increase of two orders of 
magnitude in less than five years. In addition to affecting the data 
rates, increased sensor resolution affects the following areas: 

o Greater processing requirements for image correction, formatting, 
and information extraction. 

o Increased navagation accuracy requirements to provide the necessary 
image distortion coefficients. 

o Increased archiving requirements to accept the larger data volume. 

Increased resolution of the science sensor is a positive trend for ap- 
plications. However, the acquisition of these data is totally nonde- 
terministic in that nothing is known about the quality, content, or lo- 
cation of the imagery prior to or even shortly after it is obtained. 

As a result very little' of the information is used by the scientific 
community because of such undesirable effects as cloud coverage, un- 
wanted scene content, or exposure dates that do not coincide with those 
desired. In fact, less than IZ cf all previously acquired remotely 
sensed data have been examined b> the user. 

In addition to the data deluge problem, it is currently not feasible to 
exploit Landsat data for such real-time applications as forest fire de- 
tection and monitoring and flood detection. In fact, with the excep- 
tion of applications in which the observables do not change dramatical- 
ly with time (such as oil exploration), Landsat 's usefulness is limited 
from an operational standpoint because the data are already stale by 
the time the user gets them. 

User requirements are another important trend associated with remote 
sensing missions because in the long run these requirements play a 
major role in defining the mission. Current data dissemination tech- 
niques are limited to a several-month turnaround. Many users require 
this time to be cut to hours, and a set of future users are requesting 
real-time control of the science instrument. The primary limiting fac- 
tor in data turnaround time is the processing required for image regis- 
tration. It is important to emphasize that not all users require 
real-time instrument control but the trend indicates a strong require- 
ment for faster turnaround. For applications in which the primary ob- 
servable does not change dramatically with time, or when the extraction 
of meaningful information is not time-critical, automated archiving 
techniques may be used in place of real-time instrument control. 

^References are given in Chapter VII. 
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The final trend that must be overcome to realize future missions is the 
rapidly escalating cost of ground support* The primary drivers of this 
trend are manpower costs and sophisticated equipment* A detailed ex- 
amination of the problem reveals that the following areas are the most 
critical limitations of current remote sensing missions* 

o Remote sensing missions do not employ any automated techniques of 
data evaluation for either simple annotation or data negation* A 
simple determination of percentage cloud coverage in a scene could 
effectively eliminate on the order of 50% of the data acquired* 

o Users are subject to the schedule of the mission rather than the 
mission being tailored to the specific acquisition requirements of 
the user* This results in the acquisition of tremendous volumes of 
unwanted imagery either due to poor quality caused by atmospheric 
effects, the geographic location, or the time of acquisition* 

o All data, including tracking data from remote stations, must be 
sent to a central facility for image correction, annotation, and 
packetization, and then sent to a second facility for archiving and 
eventual dissemination to the users* This results in an unneces- 
sary data link and awkward procedures to ensure critical timing 
requirements are met* This process is neither cost nor time 
effective. 

o No advantage is being taken of the onboard navigation capability to 
limit the magnitude of image distortions* By providing onboard 
control, the magnitude of the image distortions and hence the mag- 
nitude of the processing involved can be dramatically reduced. 

Automation promises to be a key driver in overcoming the many limita- 
tions associated with current remote sensing missions* Areas that 
could benefit from this technology include: 

1) Spacecraft control; 

2) Spacecraft comand generation and verification; 

3) Data unpacking; 

4) Image correction; 

5) Data archiving; 

6) Data dissemination* 

The most critical issue associated with automation technology seems to 
be whether the data processing should be performed aboard the space- 
craft or by ground support* The answer to this issue is complex and 
will involve many tradeoffs, but in the end it boils down to a question 
of economics and the long-term benefits and limitations of each 
approach* 
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In summary, the two most critical limitations of existing remote sen- 
sing missions are the deterministic acquisition of data, and the rapid 
correction of the data for radiometric effects and image distortions. 
Ihri purpose of this contract was to investigate a new approach to re- 
mote sensing that would meet future mission requirements by overcoming 
these problems. 

The primary sources of image distortion are sensor-peculiar errors, 
viewing perspective, coordinate transformation errors, and space- 
craft-induced errors. With the development of the multilinear array, 
the primary sensor-caused distortions will be the individual placement 
of the detector elements, optical effects, and orientation of the array 
relative to the sensor coordinate frame prior to flight. These errors 
remain fairly constant over time so the. resulting distortions are de- 
terministic. 

In referring to Figure I-l, viewing perspective is a well-known 
function of local earth radius and vehicle altitude and can be computed 
as shown. It is assumed that the local earth radius is known so this 
distortion is also deterministic. 

The primary error source remaining, therefore, is spacecraf t-inducea. 
The spacecraft error sources can be categorised as: 

1) Attitude determination, 

a) Star tracker accuracy, 

b) Star tracker configuration, 

c) Knowledge of star tracker misalignment, 

d) Error in star catalog, 

e) Gyro noise, 

f) Uncertainty in gyro bias, nonorthoganality and misalignment, 

g) Numerical accuracy; 

2} Ephemeris prediction, 

a) Global positioning system (GPS) accuracy, 

b) Dynamic model accuracy, 

c) GPS update interval, 

d) Numerical accuracy; 

3) Transformation error between inertial and earth-fixed coordinates, 

a) Knowlege of UTl, 

b) Knowlege of earth preccession, nutation, polar wander, and 
tidal deformation (note these error sources are generally 
included in GPS accuracy), 

c) Numerical accuracy; 

4) Misalignment between sensor and body coordinates, 

a) Knowledge of linear array orientation, 

b) Accuracy of thermal distortion model, 

c) ) Vibration modes between two coordinates, 

d) Calibraton technique and frequency, 

e) Numerical accuracy. 
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= 18O-(0 +<}>)» sin <l> 

R 

c 

== surface distance of element from the center of the field of view 



= effective resolution of i^^ element * 1 

distortion of i^^ pixel « 

where R ^ is the desired resolution, 
ref 

Fiijux^e 1^1 Viewiyuj Pex^aveative Dietortione 



For tbo i«k« of dlocuooion* asouM that all tha arror ia dua alaply to 
tha attitude dataraination ayataa* To aehiava tba temporal ragiatra- 
tion raquirament of 15 matara, it mill be nacaaaary to predict attitude 
to witbin 4 arc-aaconda aa illuatratad in Figure 1~2. Tba accuracy of 
tba atata-of-tba~art syataM uaing tba HA8A atandard atar tracker and 
gyro ia IS arc-aeconda (2o) (Kaf 2). 



Figure 1-2 

Error Budget for RegietraHon Aoouraay of IS m 


Evan with tba davalopmant of charga-coupled davica (CCD) star trackara« 
tba attituda datarmination capability will be around 6 arc-seconds 
(2o)* Note that 2anuabera have been used bera corresponding to 95X of 
tba data. If lo numbers corresponding to 67X of tbe data are uaed» tbe 
accuracy goals can be met. However, by adding even one more error term 
sucb as misalignment between senaor and body coordinates of 2 arc-sec- 
onds (2-axis alignment accuracy acbievable with optical alignment 
cubes), tbe Cotel error budget is exceeded. From tbe previous discus- 
sion, which ignored many error sourcea* it ia clear that another ap- 
proach is required. 

Tbe key to real-time image correction lies in being able to accurately 
determine tbe location of the science sensor's boresight in earth-fixed 
coordinates. Witb this capability it is not only possible to provide 
real-tiiBS aieasursaents of tbe image distortions, but with the advent of 
tbe multilinear array (MLA) it sMy be possible to provide real-time 
image correction using either resampling or special design of tbe MLA 
focal plane. This capability also provides tbs heart of a pointing 
system capable of deterministically acquiring isiagery at specific 
eartb -fixed coordinates. 

Shortly after definition of tbe feature identification and location 
experiment (FILE), Martin Marietta began tbe development of s landmark 
tracker or GCP detector centered around experience gained witb terminal 
guidance systems. Tbe primary purpose of tbe landmark trcsker is to 
provide periodic awasurements of the science sensor's boresight posi- 
tion to be used as input to a navigation system. Previous studies (Ref 
4 thru 7) have sbotm that tbe landatark tracker cannot solve for both 
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position and attitude without aupplemental neaaurementa from another 
source. Another limitation of the landmark tracker operating in the 
visible spectrum is that observations can be obscured by clouds and the 
correlator will lock onto a false target. For these reasons, the con- 
ceptual navigation system consists of the landmark tracker, a CPS re- 
ceiver, and two NASA standard star trackers to replace attitude meas- 
urements when the landmark is obscured. 

The purpose of this contract was to simulate operation of a navigation 
system centered around image correction, and analyse performance of the 
system under a variety of conditions. Also of interest was a sensi- 
tivity analysis to determine the optimal sampling interval for each 
sensor, and the total number of landmarks required to satisfy mission 
requirements. 

In the past there has been a lack of coordination between the scien- 
tific user community and the engineers responsible for spacecraft de- 
sign. This has resulted in a physical separation between the design 
and implementation of the science payload and the control system. An 
example of this thinking is shown in the multimission spacecraft (HMS) 
where subsystems are treated as modules and the payload itself is 
physically separated from the control portion of the vehicle. This 
type of design, although attractive from a standardisation point of 
view, ignores the inherent relationships between user requirements, 
spscecraft control requirements, and ground support requirements. It 
is possible that rather than saving significant costs, standardization 
may result in higher end-to-end costs. This design philosophy must be 
reevaluated with regard to future missions. In the end-to-end design 
of remote sensing spacecraft, the primary emphasis in the guidance and 
control system must shift from simply estimating the ephemeris and at- 
titude of the spacecraft to estimating the position of the science sen- 
sor's FOV on the earth's surface. This provides the basis for both 
real-time image correction and deterministic data acquisition through 
sensor pointing. This shift of emphasis will impact the design of the 
entire spacecraft. For example, if the science sensor boresight posi- 
tion is to be determined, it is desirable to place the gyro package 
close to that sensor to reduce the misalignment between the two. The 
(MHS) configuration, which provides a physical separation between the 
payload and the guidance and control system, may not satisfy the re- 
quirements of many future remote sensing missions. 

This report has been separated into four volumes. The first volume 
entitled "Executive Summary" is intended to provide an overview of the 
iisage correction system, a brief examination of the analysis, and a 
discussion of the impacts on future missions. This volume, entitled 
"Analysis and Simulation Results," provides a detailed account of all 
work performed under the contract. Volumes 111 and IV contain a 
detailed description of all software developed and used in this 
contract. 
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This volumt has bean separated into the following chapters: 


I. Introduction; 

II. Onboard Isuige Correction System Overview; 

III. CCP Navigation System Simulation; 

IV. CCP Analysis; 

V. Image Correction Accuracy Measurement; 

VI. Conclusions end Hecosssendations; 

VII. References* 
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ONBOARD IMAGB OORRBCTXON SYSTEM OVBRVUM 



Tb« oaboArd iMg« eomeeioii tyctm has baan aapaMtad into thraa pri* 
■ary aubayataaa aa illuatratad in Figum XX>1* Tba navigation procaa~ 
aor aecapta aaaauraaanta froa tba GPS raealvart tba atar t<aekara» and 
tba landaark traekar, and aolvaa for spaeaeraft atata* Tba rogiatra** 
tion procaaaor ia rtaponaibla for eorralating liva iaagary nitb atorad 
rafaranea ground control points to produea a landaark sighting vaetor. 
Ibis procaaaor is also raaponaibla for providing ■aaauroaanta of tba 
distortion coafficianta« and for providing iaaga raaaapling if tbia 
oparation ia parfonaad on board* Tba pointing nount controllar baa 
diffarant uaaa dapanding on tba application. If raal-tiaa iaaga cor> 
raction is baing providad through tba uaa of a curved focal plana nul* 
tilinaar array (KLA)i tba pointing a^unt would ba uaad to anaura that 
tba inatruaant borasigbt always pointed at nadir* It aalactiva acqui<> 
aition ia baing provided, tba pointing aount could ba used to provide 
additional flexibility* Tbia chapter daacribea each of tbasa aubaya- 
teas in graatar detail* 



Hgux'e IJ-1 Onboard Image Correction Syetm Block Diagram 
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A. REGISTRATION PROCESSOR 


Figure 11-2 illuetretee the functional flow of the registration proces- 
sor. A digital representation of all the ground control points is 
stored on a high-density tape recorder (or bubble steaory) and read into 
randoa access Mnoty as required. Because of the reliability consider- 
ations associated with the recorder, a nusd>er of GCPs will be read into 
RAN at one tine. As the spacecraft proceeds in orbit, the landmark 
that is to be viewed next is determined and the digital reference image 
is extracted from RAN. Information stored in the GCP file, along with 

uncertainty in vehicle position and attitude, helps to define the loca- 

tion and sise of the search area. For an operational system it is ad- 
visable to use a fixed-sise search area. Since a star tracker is being 
used for backup, the sMximum attitude error will be in the vicinity of 
IS arc-seconds (2 o ) and the global positioning system will produce an 
error around 12 meters. Therefore a search area that is 32 pixels 
larger than the GCP sise will provide far more area than required. For 
this reason, and the fact that a GCP sise of 32 produces Che best re- 
sults, a search area sise of 64 was chosen. 

Proceeding with registration processing, a test is initiated to deter- 
mine if the location of the present scan line is coincident with the 

position of the search area. If it is not, a new satellite state is 

obtained from the navigation proceasor. Ibis new state is then trans- 
forsied into sensor pointing information for use by the pointing mount 
controller for several purposes. First, if a curved focal plane array 
has been used to provide real-time image correction, the pointing mount 
suist isolate the array from the spacecraft isotion to ensure that the 
central element of the array is viewing directly along nadir. If the 
pointing mount is being used to provide deterministic data acquisition 
according to earth-fixed coordinates, the boresigbt position provided 
by the registration processor will be used as a reference tu point 
from. Following the determination of boresigbt position, the loop con- 
tinues until the current scan line is coincident with the search area. 

At this point, image scan lines extracted from the sensor are stored in 
a buffer and become the search area data. The search area is then 
registered with respect to the GCP and the registration vector passed 
to the navigation processor, which uses this information to update the 
estimate of vehicle state. The entire process is then repeatei for the 
next landmark. The following paragraphs describe in more detail the 
operations performed during each of these functions. 

1. Organisation of the GCP Pate Bese 

The architecture of a GCP data base is closely related to the require- 
isents placed on a particular mission. For the advanced Landsat-type 
mission sMdel being considered, we selected a data base that is primar- 
ily a sequential access file but allows impleisentation of random access 
for missions requiring deterministic data acquisition through sensor 
pointing. Ibe data base has been broken into four levels of directory 
as shown in Figure II-3. 
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Hgure II - Functional FUm of Regietration Prooeaaor 

The first level contains one entry for each orbit of the satellite. 

For Landsat, this would nean 16 days worth of entries since the orbit 
is repetitive with this period. Each entry in the orbit table can be 
broken down into sections of an orbit where the size of the section 
will be determined by equating the number of GCPs contained in a sec- 
tion with block size to be accessed from the recorder. Finally, an en- 
try in the orbit section table can be broken down into the GCPs con- 
tained in that section. This final level of data base is called the 
count table and contains all the necessary information to locate the 
position of the search area plus the starting address in RAM containing 
the digital representation of the GCP. 


Orbit No. Orbit Section Count Table 



Landmark 1 

Start 
Counter 1 

Stop 

Counter 1 

Element 
Start 1 

Element 
Stop 1 

Address 1 

Landmark 2 

Start 
Counter 2 

Scop 

Counter 2 

Elenant 
Start 2 

Element 
Stop 2 

Address 2 







Landmark N 

Start 
Counter N 

Stop 

Counter N 

Element 
Start N 

Element 
Stop N 

Address N 


Landmark Data File 
Latitude 

Longitude 

Landmark Site 
Mean Radiance 
Sun Angle 

Correlation Convergence 
Indicator 

Future Flag Words 

Video 

Data 

EOF 


FiijuPo Organization of the GCP Data Baae 
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Determination of Sensor Boreeight 


The position of the science instrument boresight in earth-fixed coor- 
dinates is required for pointing control (Fig. 11-4). In referring to 
Figure II-4, the boresight position may be found by forming a unit vec- 
tor along the boresight in the sensor coordinate system 

Fol 



wiiere is a unit vector in landmark tracker coordinates. This 
vector may then be transformed into inertial coordinates through 

i!l ■ iTg bTl Ul 
where 

JUl ■ unit vector along Che sensor boresight in inertial coordinates, 

gT^ " transformation from sensor coordinates to body coordinates 

accounting for pointing angle offsets, 

iTg ■ the direction cosine matrix transformation describing the 

attitude of the vehicle* 
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Since 

R - P - U 

Nbere 

R • radius vector from the origin of tbe earth -centered inertial 
coordinate system to tbe boresigbt surface position^ 

P * space craft position vector in ECI coordinates* 

M, " measurement vector from tbe spacecraft to tbe boresigbt posi- 
tion on the earth's surface* 

It is necessary therefore to define ll* This may be done using tbe laws 
of sines and cosines as follows* With |R|, tiie magnitude of R known* <t> 
may be determined by tbe following process* ~ 

J&J liL. yd 

sin e sin C sin (180-C) 

Where 6 is determined by 
U*P 

® IHTF' 

Solving for C in terms of 6* R and ^ produces 

Ip| sin e 

C = 180* - sin-1 • 


Since 

+ c + e = 180** 


|p| sin e 

4. « 180“ - C - 0 = sin-1 ^ 


Using the law of cosines 

1m| = (|mP)^ = (|rP + llP - 21 p1|r| cos 

Ip| sin e , 

= (|r|2 + |p|2 - 2|p||r| cos (sin-1 — ^ e))\ 

Since M is directed along U* 

M = |m1 U 

and 

R - P - M* 

The longitude (L) and latitude (A) of the surface position may be de- 
termined if R is transformed into earth-fixed coordinates. This trans- 
formation, gTi, is a function of time, earth rotation, polar wan- 
der, nutation, etc and is assumed to be known to the system. 
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Therefore 


■ E^I St • 



Ri ■ the projection of Re along the i axia» for i * x» y» s and 

The respective earth-fixed coordinate axea conaiat of a along the 
earth's nominal spin axis, x from the earth 'a center through the inter 
section of the equator and the Greenwich meridian, and T completea the 
right-band set. 

The longitude and latitude of the boresight aurface projection are 
given Iqf 


A = sin 


-1 - 5 - 


L » sin" 


R cos A 


Extraction of Search Area Data 


The location and sise of the search area is coa^letely defined by the 
line start, line stop, element start, and element stop inclcators con- 
tained in the count table as shown in Figure II-S. 

Present Sensor 
Scan Line 


Field of View Over Time 


Element 

Stop 


Start Count (i) 




Clement 

Start 


Search Area (i) 


Fijux^ II-5 Explanation of Papametere in Count Table 
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The two long horizontal lines enclose a diagram of the imaging sensor's 
field of view over time* The vertical swaths depict one scan of the 
sensor's instantaneous field of view, and the small squares within the 
scan line are individual picture eleswnts. The start count contains 
the number of scans between the present scan line and the first line of 
the search area. The strp count contains the size of the search area 
in scan lines. Once it is known that the search area is within the 
field of view, it is necessary to determine its location. The element 
start count contains the distance, in picture elements, of the search 
area from the start of the scan line. The eleoMnt stop count contains 
the distance to the end of the search area from the start of the scan 
line. These four parameters totally describe the position and size of 
the search area relative to the current position of the sensor's field 
of view. For a constant search area size, two of the parameters may be 
eliminated to reduce the onboard storage requirements. The center of 
the search area is coincident with the estimated center of the landmark 
and its size is a direct function of the uncertainty in satellite posi- 
tion and attitude* To accommodate the varying size of the search area 
over time, the parameters in the count table can be changed. 

The algorithm implemented to extract the search area from the sensor's 
field of view (FOV) is shotm in Figure 11-6. The start count indicator 
is decremented each time a scan line is read by the sensor. The decre- 
ment process is driven by the same clock signal that drives the scan 
mechanism so the counter does not diverge. When the counter reaches 
zero, the search area has been encountered and the second loop in the 
figure is entered. In this loop image data from the sensor are read 
into a buffer each time a scan is initiated. The buffer does not con- 
tain the entire scan line, but only the data between the element start 
and element stop flags previously described. After a scan line has 
been recorded, the element stop count is decremented. When this count- 
er is zero, the entire search area has been read and the correlation 
pvocess can begin. 

4. Registration of Sensor Data 


The registration processor is centered around a sequential similarity 
detection algorithm (SSDA) first identified by Barnea and Silverman 
(Ref 8). Various correlation algorithms were analyzed under a separate 
contract (Ref 9) and the results of this study indicate that the SSDA 
can be used to detect features to within one-tenth the resolution of 
the science instrument. In addition the SSDA was found to be much sim- 
pler to implement in an operational system primarily because it con- 
sists of simple additions whereas other algorithms involve multiplica- 
tions and powers of large numbers. 

A ground control point may be defined as a small area within the sen- 
sor's FOV whose characteristics ate easily detected using a registraton 
algorithm such as the SSDA. It is not necessary to process imagery 
from the entire FOV, but only an area whose size ensures that the GCP 
will be located within f.t« boundaries. The size of this search area is 
a direct function of the uncertainty in the vehicle's position and at- 
titude. However, for an operational system where timing requirements 
are critical, it is adviiiai 'e to determine the maximum expected error 
and to fix the search area sizt'. 
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Fujux'e II-6 Seai'ijh /Ij't-’U AUjai'Wm 

Tbe registration process SMy be described in tbe following way* Let 
the search area be defined as an LxL area of digital picture elements* 
Tbe isuigery nmy then be defined by a function* 8, that describes the 
gray scale* or recorded radiance* in relation to position coordinates* 
i*e* * 

8(i*j) - 

vAiere is the gray scale of the i,j^^ picture element of the 

search area defined by 

l<(i*j)<L. 

Let the ground control point be defined similarly aa an MxM area with 
an image function 

G(l*m) - 
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where Ri _ is che grey scale of the picture element of the 

GCP defined by 

l<(l,m)<M. 

A subimage (Fig. 11-7} of the SA may be defined as an MxM area whose 
upper left coordinates (n,o) lie in the range 

l£<n,o)£L - M + 1. 


X I GCP 



Figure II-? Subimge Definition 


A subimage whose upper left coordinates are n,o will be referred to as 
the n,o^b reference point. 
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Th< s«nsor data is ragistarad by maaauring tba sijsilarity batwaan aaeb 
MXM subiaaga vitbin tba saareb araa and a digital rapraaantation of tba 
GCP storad onboard. Iba rafaranca point that produeaa tba bigbaat da~ 
graa of aiaiilarity with tba GCP ia tba bast ragistration of tba GCP 
and can ba labalad with tba aama aartb fixad eoordinataa as tba GCP. 


Tba 88DA algoritbn nay ba iaplanantad to dataet aiailarity batwaan a 
rafaranca point and tba GCP tbrougb tba following aquation: 

Similarity * 

Where 

|no* n»otb iteferanca Point, 

Mean value of tba aubiouige located at tba no^b referanca point, 

G • Mean value of tba GCP 

Ibe entire registration process can then be described by tbe algoritbai 
shown in Figure 11-8. 


^ Trial Poaitlono 



Figure II~8 

Registration of Sensor Data via GCP Detection 


Sooie evidence suggests tba naan value of separata features should ba 
independently computed. For axampla, tba mean radiance of vegetation 
wnitld be computed separately from bare earth, and this naan tiould only 
be used to norsialize vegetation pixels witbin tbe subimage. 


Several approaches to subpixal ragistraton ware evaluated under inde- 
pendent rasearcb. Iba first approach consisted of interpolating tba 
imagery corresponding to tba referanca point that produced tba highest 
degree of similarity. Tba interpolation of tbe isMgary consisted of 
resampling tba imagery at a 10:1 scale. Tba resampled image was than 
registered using tbe 880A algorithm with tba GCP. Ibis technique 
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proved to be useful in providing a registration accuracy to within 0*1 
Pixel (Ref 10). However, Che processing involved in resampling Che 
image was found Co be excessive for an onboard syscem. Therefore a 
second approach was considered. 

As sisiilaricy is being measured by Che S8DA Che correlacion surface is 
scored in a separace array. When Che besc fic is located, Che surface 
is inCerpolaCed by ficcing a cubic spline, and Che minimum value of Che 
funccion is obcained. Resulcs obcained show chac Chis Cechnique is 
also capable of providing a regiscracion veccor with an accuracy of 0.1 
Pixel. The processing involved in interpolating Che surface is also 
well within the capabilities of an onboard syscem. 

Cloud Supression 

A small percentage (lOX) of cloud coverage within the search area can 
affect the correlation process in such a way Chat the GCP is no longer 
recognisable. This effect can be minimised by simply eliminating every 
picture element representing a cloud from Che correlacion and normalis- 
ing the coefficient according Co Che number of pixels correlated. Pre- 
liminary resulcs indicate Chat 20 to 40X cloud coverage can be toler- 
ated depending on the type of landmark and what portion of Che scene is 
obscured. The problem then becomes one of detecting clouds in the 
scene onboard in real time. 

The feature identification and location experiment (FILE) system pro- 
vides the capability of classifying the sensor field of view into the 
four generic features of bare land, vegetation, water and clouds. Fig- 
ure II-9 illustrates how various feature types can be separated using 
two sensors operating at 0.63 and 0.83 pm. 

Although the polygons in the figure overlap to a certain extent because 
they represent 99X confidence limits, the majority of cases can be ex- 
pected to fall in unambiguous points on the plane. This technique was 
incorporated into the concept for onboard GCP detection to provide a 
pixel-by-pixel suppression of clouds. 

False Lock Detection 


Since the registration vector is being used to update the vehicle navi- 
gation filter, it is of the utmost importance to detect when the corre- 
lator has locked onto a false target. Analysis indicated that correla- 
tion surfaces resulting from false target lock always converged more 
slowly chan Chose resulting from a correct lock. The first approach to 
false lock detection was to measure the convergence rate of the surfa.*e 
and implement a threshold discrimination algorithm. However, the con- 
vergence rate was found to be a complex function of several variables 
that included target contrast and geometry, indicating that a simple 
threshold algorithm would not perform well for all cases. A more thor- 
ough analysis showed that a better indication of false lock would be to 
measure the similarity between the correlacion surface obtained from 
registration of the search area with respect to the GCP and the corre- 
lation surface obtained with GCP autocorrelation. The false lock coef- 
ficient is measured using a classical correlation algorithm. 
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Figure II-9 99% ConJ'idence PoUjgone for Feature Claeeification 

The accuracy of false lock detection was studied under the analysis and 
results indicate that 99X of false locks may be avoided. 

Correlator Design 

In an operational system, it is highly desirable to implement much of 
the registration processor functions in a hardwired system. This sec- 
tion discusses one possible implementation of the correlator. 

There are four elements in the landmark correlator: 

1) Data selection; 

2) Data buffering; 

3) Correlator; 

4) System control. 


The following paragraphs describe each of these elements. 







• • P«t« 8<l»ctioB - Dtu ••Ueelon is • two-tt«p pip«lliM op«r«tieo* 
Data ara aaauaad to ba g*Mrata4 by tba aanaor aa aa 8~bit wida paral- 
lal data atraaa* lacb 8 bita dafina tba valua of ona pixal* t¥o »yn- 
cbronisatioa aignala ara alao aaaua«d to ba available from tba a«i\aor* 
Ibaaa ara OidTA CLOCK and baglnniag of lino (BOL)* DATA CLOCK la a 
pulaa that indieataa that tba data linaa contain data* BOL indicatas 
tba tiaa at vbicb tba aanaor ia atarting a awaap or acan* 


Pigura 11-10 illuatrataa tba paraaatara involved in tba data aalaetion 
procaaa* Tba area in vbieb tba landaMrk raaidaa ia amallar than tba 
aanaor field of view* Tba landmark area ia defined LIKE CODIIT by PIXEL 
OOOMT, or nuabar of linaa by nuAbar of pixala. Kacb line ia aaauaad to 
contain tba aane nuad»ar of pixala. Tba ayabolic upper left-hand comer 
of tba landmark area ia defined by tba valuaa LIMB BIAS by PIXEL BIAS. 
Some aynebroaisation point auat be aupplied for tba LIME BIAS maaaure- 
mant. Tbia aynebronixation point could be aa aimple aa a timer aup- 
pliad by tba landaMrk ayatam control element. Tba aynebroniaation 
point for PIXEL Bl\S ia tba BOL aignal generated bP tba aanaor* 


BOL tensor Field of View 
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Fijui'e II~10 Patvaneteva Involved in Data Selection Pvoaeaa 

Tbe firat atep in tbe data aelection procaaa ia to aelect tbe proper 
acan linaa. Thia ia accompliabed by initialiaing two countera with tba 
valuea of LIMB BIAS and LIMB COUNT. Bacb time a BOL aignal ia gener- 
ated by tba aanaor, tbe LIMB BIAS counter ia decraaMinted by one. Nben 
tba LINE BIAS counter raacbea aero, tbe aignal LSTKT ia generated, in- 
dicating tba firat line of tba landawrk area* At tbia time tba counter 
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eonttining LIMB COUMT will b« dwerMMnttd for ooeb BOL pulto* Wboa tbo 
LIMB COUNT roaebot Mro, L8T0P la ganaratad. A aiapla laCeb la util- 
itad to control an L 8 AVB aignal that danotaa tba tiaa fron LSTK to 

L8T0P. This ia tba tina aban aanaor acan linaa corraapond to tba land- 
aark area* 

Iba LSAVB aignal anablaa tba aacond atap of tba data aalaction proc- 
aaa. Tbia atap worka in tba aaM way aa tba firat. Inataad of count- 
ing linaa, bowavor, pixala from BOL ara countad. Tba output of tba 
5 alact proper pixala atap la tba PBAVB aignal indicating that pixala 
vitbin tba landmark araa ara being ganaratad by tba aanaor. Tbia PBAVB 
aignal ia uaad to drive tba next landmark correlator alamant, tba data 
buffering element * 

P*ta Buffering - Aa would be expected, data buffering involvaa both 
an input aLaga and an output ataga. Tba input ataga ia atartad by tba 
PBAVB aignal. Mtan tbia aignal ia true, da ta from tba aanaor ara 
ab if tad into a BAM buffer in a direct memory accaaa (DMA) faabion. Tba 
DMA logic ia implamentad by two down-countara, one counter controlling 

tba RAM tranafar addreaa and tba other counter controlling tba tranafar 
count . 

The output ataga of data buffering nay begin before the entire landmark 
area data have bean collactad bacauaa tba correlator alamant will oper- 
ate in a aerial pipeline manner. Unfortunately tba order in which daU 
are input to the RAM buffer ia not tba aama order in wbicb data ara 
output from the RAM buffer. Por example if tba landmark araa ia diman- 
aioned aa 40 linea by 40 pixala, tba incoming data atraan ia 

LlPl. L 1 P 2 , L 1 P 3 , ..., L 1 P 39 , LiP 4 o, L 2 Pi, L 2 P 2 L4oP40‘ 

If the landmark ia 16 linaa by 16 pixala, tba data atraam output from 
the RAM buffer will be 

L1P2, ..., LiPjs, L^Pig, L 2 Pi» L2P2 LigPjg 

CO compute Mb (1:16, 1:16) and C(l:16, 1:16), and 
L1P2, L1P3, ..., LjPjg, LiPi 7, I*2P2» •••» LigPi 7 

to compute Wm (1:16, 2:17) and C (1:16, 2:17), etc. Ueb output data 
atraam from tba RAM buffer may be bandied in a DMA moda but tba ayatam 
control element muat continually meat tba initialiaation logic of tba 
output ataga to anaura correct data aalaction. 

c. Correlator - Tba correlator ia implamentad aa a four-atap pioa- 
linar %ie required atapa ara: ^ 

P P w 

Step 1: Compute (3S ■ — ii» 

n - m 


Step 2: Co^Hite Cx (i,j) • Vm * Lk (i,j) 

where Lx (i,j) « La (i,J) - lai. 
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Step 3: Coaputc Of (l»i) * Mb (k,^) “ Cx (i*J)i 

8t«p 4: Compute C ■ C ^ Ppl(l»J)l* 

Vitb ttoii opprooeb it is noeoooofjr tbot otop 1 k« eomplotoly finUbed 
beforo tbo otbor itope eon bogio* Oneo tbit it <ioo«« boutvor* ttopt 2» 
3, and 4 ouqr oparata tiattltanaeutly witb atap 3 oaa data point babind 
ttap 2 and atap 4 ona data point babind atap 3t Tba output of atap 4 
it tba eorralation aalua for ona poaition of tba iandnark witbin tba 
landmark araa. Computation of all tba eorralation valuaa ia tiad in 
vitb tba output ataga of tba data buffering alamant and tbia ia da- 
acribad more fully in tba following paragrapb. 

d. Syatam Control - Baeauaa of parfoimanea raquiramanta tba tbraa ala- 
BMnta juat ^aacrikad (data aalection. data buffering and tba correla- 
tor) abould be implamentad uaing bardwirad logic. Hardwired logic 
could alao be uaad for tba ayatam control alamant although tbia appaara 
neither neeeaaary nor daairabla at tbia time. A >^ucb more attractira 
approach could be tba uaa of a dedicated mieroeoaqnitar. Tba ayatam 
control functiona can be daacribad quite aaaily in aoftwara and indiea- 
tiona are that obtaining aufficiant parformanca would not be a prob- 
lem. The ayatam control functiona required are: 

1) Data buffering input ataga initialisation; 

2) Data buffering output ataga initialisation; 

3) Correlator atap 1 and 2 synchronisation; 

4) Correlator step 4 data acceptance. 

As data is collected from atap 4 of tba correlator tba ayatam control 
alamant must keep track of tba minimal C value and where tbia value 
occurs. This information ia than fumiabad to other procaaaaa for 
datarmining spacecraft state. 

a. Landmark Correlator, karformanca Hatimata - Figure 11-11 abowa a 
top-laval block diagram of tba alamanta Just~diacussad. Several ap- 
proaebes ware considered and aoaw ware iaadiataly discarded because 
sufficient performance could not be obtained. Other approaches involv- 
ing more parallel operations could be used to obtain even higher per- 
formance but cbe additional bardwara complexity did not appear justi- 
fied at tbia time. 
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-11 Correlator Design 









Two types of tecbnology wore considered for berdwere implenente~ 

Cion — TTL end CMOS/SOS. In botb coses Cbe liaiCing perforaance eleaenC 
wes step 3 of tbe correlator* For tbe TTL iapleaentation a a a ni aua 
clock* of 10.5 MHz could be obtained. If seaicustoa GM0S/S08 LSI is 
used, this seas eleaent has a siaxiaua clock* frequency of 12 MBs and 
has ttie additional advantage tbat it requires significantly less power 
and coaponent count. 

With tbe step 3 correlator liaited to 12 MBs, one correlation coeffi- 
cient can be coaputed every 22 us for a 16x16 landaark. If Cbe land- 
aark area is 40x40, correlation of tbe entire area requres 11 us. 

Since tbe original goal was to do tbe whole correlation in less than 1 
second, this perfonsance level was considered adequate* 


NAVIG4TI0B PROCESSOR 


Figure 11-12 illustrates the functional processing performed by the 
navigation processor. The most critical function performed is the 
sequencing of measurements and supervisory control of the registration 
processor. The sequencer first decides, based on a measurement pro- 
file, what type of measurement should be made and when the measurement 
should be made. The types of measurements the sequence tables estab- 
lish are star trackers 1 and 2 and GPS. After determining the time of 
the next measurement, the sequencer decides, based on Che total number 
of scan lines left Co the landmark and Che scan rate of Che sensor, 
whether there is enough time before the landmark sighting. If where 
is, Che measurement is processed normally. If there is not. control is 
defaulted to Che registration processor until Che registration vector 
is returned. 

After ti»e aeasurenent is obtained, tbe estimated state vector is inte- 
grated forward to Cbe measurement time. Ibis process is part of tbe 
extended square root navigation filter but is implemented in a separate 
package to siaq>lify control* Ibe raw measurements are processed to 
compensate for knowlege of sucb error terms as bias and misalignment. 
These compensated measuresmnts are used as input to tbe Kalman filter 
package tbat then estiisates tbe spacecraft state. Tbe spacecraft state 
includes: 

1) Position; 

2) Velocity: 

3) Attitude; 

4) Gyro Drift, nonortboganality, and scale factor; 

5) Science sensor misalignment. 


*Maximum clock frequency derated 50Z for spacebome environment* 
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Pigux*e II-12 Navigation Funotional Deeanption 

The navigation procesaor is raaponaibla for aending tba navigation 
state to the registration procesaor at the beginning of each scan 
line. This information ia used to compute the distortion coefficients 
and to control the senaor pointing mount. Attitude is propagated be- 
tween measurements by a HASA standard gyro package. The pulse train 
output from the gyros is compensated for knowledge of bias, nonortbo- 
ganality, and scale factor and then integrated to provide vehicle 
attitude. 

Bach of these processes is described in snre detail in Chapter IV under 
the sections discussing dynamics and integrator, sequencer. Kalman fil- 
ter, and gyro processing. 
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C. POINTING MOUNT CONTROLLER AND LINEAR ARRAY 
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The pointing smunt controller can be implemented in two different ways, 
depending on the configuration of the image correction system. If a 
curved focal plane array is used, the pointing mount controller would 
be responsible for maintaining tbe boresigbt of the instrument along 
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nadir« In a system where image correction is provided through resampl-* 
ing» the pointing mount controller could be used to provide determinis* 
tic data acquisition by pointing to specified earth**£ixed coordinates* 

In a curved focal plane configuration, the pointing mount controller 
would form an error vector that represents the angular offset between 
the boresight of the linear array and nadir* This error vector would 
then be used as an input to the pointing mount to correct for the 
error. In this way, the viewing angle of the MLA would be isolated 
from the spacecraft motion caused by either attitude limit cycles, or 
attitude corrections* 

In a selective acquisition mode the sensor boresight position would be 
compared to the desired pointing angle to generate an offset that would 
be used as an input to the pointing mount. Again the viewing angle of 
the sensor can be isolated from the spacecraft motion to simplify the 
image correction process* 

The idea behind a curved focal plane array originates from the under- 
standing that, with the development of the MLA, the primary image dis- 
tortions will be caused by the viewing angle and the curvature of the 
earth (Fig. II-l)* This nonuniform sampling is relatively constant and 
is symmetrical' about a vector pointing along nadir* It was therefore 
felt that constructing an array with nonuniform sampling for each ele- 
ment that compensated for these effects could eliminate along-scan re- 
sampling, which is one of the largest bottleneck of remote sensing 
(Fig. 11-13). 

The limitations of this approach are that the central element must al- 
ways point at nadir to within 0.5 the resolution angle of the sensor. 
While this is anticipated as being feasible with either Gimbalflex or 
the angular suspension pointing system, a constraint is placed on the 
system in that it can no longer be used for deterministic data acquisi- 
tion through sensor pointing. On the other hand, the cost savings and 
increased data turnaround time make this an attractive possibility that 
warrants further investigation* The key considerations that need to be 
addressed are: 

1) Can the curvature fo the focal plane maintain a precision that en- 
sures additional distortions are not introduced? 

2) Is the variation in the local curvature of earth great enough to 
introduce errors greater than 0*5 pixel? 

3) Are variations in the orbital altitude great enough to cause 
distortions? 

4) What is the cost savings of performing real-time image correction 
using this technique versus complex resampling processors? 

3) How real is the desire to provide deterministic data acquisition 
through sensor pointing? Are the cost savings of deterministic 
data acquisition greater than the cost of producing the curved 
focal plane array? 
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Figure 11-12 Curved Foaal Plane Array 

The answers to these questions will form the basis for selecting the 
system architecture to be used on the next generation of remote sensing 
vehicles. Ho%rever, even if a curved focal plane is not chosen* the 
separation of the science instrument from the nation of the spacecraft 
through the use of a pointing mount will reduce the distortions and 
simplify the process of image correction while providing deterministic 
data acquisition. 
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111. GCP NAVIGATION SYSTEM SIMULATION 


The GCPSIM program has been configured to provide scientific simula- 
tions to predict the performance of the GCP detection system* Two 
types of simulation modes have been incorporated. The first mode sim- 
ulates the spacecraft the spacecraft environment and all measurements. 
The second mode extracts ground control point measurements directly 
from Landsat imagery. 

The program was set up to provide extensive error analysis rathor than 
simply a covariance analysis. Although covariance analysis provides a 
great deal of information, in many cases this information can be inac- 
curate or misleading. For example, there are cases where the covar- 
iance matrix converges over a period of time while the actual state es- 
timate diverges from the true state. For this reason the simulation 
models the true state of the vehicle using the dynamic equations of mo- 
tion. By propagating the true vehicle state, one can perform a covar- 
iance analysis and also generate the actual state error and measurement 
residuals. The two modes of operation are described in the following 
sections. A thorough description of the software is contained in Ap- 
pendix A. 


A. SIMULATED GCP MEASUREMENT 

GCPSIM was designed to provide the ability to analyze the effect of 
various measurement sequences. This was especially important when 
studying the effect of GCP spacing, missed GCP sightings, and the ex- 
pected accuracy after traversing a large body of water. The measure- 
ment sequencer designed for GCPSIM allows any mixture of GCP, global 
positioning system or star tracker measurements and time delays (peri- 
ods during which no measurements were made) of any length. The se- 
quencer determines the type of measurement and the time at which the 
measurement should be made. The true vehicle position state is then 
propagagated forward to this time by integrating the nonlinear equa- 
tions of motion with some additional process noise to account for mod- 
eling errors. The attitude state is propagated by looking up the body 
rates in an attitude profile table and integrating these rates. 

The true vehicle state is used along with a measurement model to gener- 
ate an ideal measurement vector. The ideal measurement is then corrup- 
ted with noise, bias, and misalignment terms and compensated for know- 
ledge of these values. This allows a careful analysis of the effect of 
misalignment on the state solution. It is important to understand the 
effect of bias and misalignment between the landmark tracker and body 
axis because this is the largest unknown factor contributing to a 
pointing error. It is possible to provide frequent inflight calibra- 
tion for these misalignment errors, but it is difficult to model, for 
any length of time, the various processes that cause the misalignment. 
For example, thermal gradients across the vehicle and vibrational modes 
within the flexible structure are complex functions of such things as 
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structural design, sun angle, physical properties of the material and 
many other factors. These processes are the most difficult and least 
understood of all engineering problems. 

The compensated measurements are used as inputs into an extended square 
root Kalman filter, which estimates the true vehicle state. The ex- 
tended filter propagates the estimated navigation state, the state 
transition matrix and the process noise array between measurements by 
integrating the various differential equations using a fourth order 
Runge Kutta Gill process. It is anticipated that a linear propagation 
would be satisfactory for S-second and CPS update intervals. However, 
since the CPS update interval was one of the parameters investigated, 
the more general approach of state propagation through integration was 
used. The estimated attitude state is propagated by a gyro model that 
corrupts the output with gyro drift, noise, nonorthogonality, scale 
factor and misalignment. The gyro output is compensated in a similar 
fashion to the measurement model. 

The estimated state is used to form an estimated measurement, which in 
turn is subtracted from the true measurement to obtain a residual. 

This measurement residual and calculated Kalman gain are used to update 
the state estimate, which can be compared with the true state to yield 
the error. The entire process continues until the spacecraft is prop- 
agated forward to the run stop time. 

The remainder of this section describes each of these processes in more 
detail, and Appendix A provides a detailed software description of the 
program. 


Measurement Models 


This section summarizes the mathematical measurement models used in the 
simulation of an onboard GCP detection system. The types of measure- 
ments used include GPS position and velocity, star tracker sightings, 
and landmark tracker sightings. The design philosophy behind the mea- 
surement models is that the actual vehicle state is used with a geome- 
try model to yield an ideal measurement vector. This ideal vector is 
then corrupted with bias, noise, and misalignment to provide the actual 
sensor output. The sensor output is then compensated for some estimate 
of the error terms. This compensated measurement then is useo by the 
filter to estimate the vehicle state. The benefit this design approach 
provides is that it enables a detailed analysis of sensitivity to mis- 
alignments and compensation abi.ity. It is also expected that the se- 
vere requirements associated with onboard image correction will require 
the onboard estimate of misalignment terms as those between the 

MSS and body corrdinates. With this approach it will not be difficult 
to modify the filter to solve for these terms. The remainder of this 
section has been separated into subsections corresponding to each mea- 
surement mode. 


II1-2 


a. Landmark Tracker Measurement Modal Development - The landmark location 
on Che surface of Che earch in local landmark coordinaCes will be a 
function of Che altitude (Aj^) above Che earth's mean radius 





0 


0 


However, in earth-fixed coordinaCes, Che landmark will have Che earth's 
mean radius (rg) added Co the altitude. Using the angular transfor- 
mation from local landmark to earth-fixed coordinates produces 


^ “ E^'L ^ " E^L 



•"V 

0 

0 



As shown in Figure III-l, Che position vector of Che spacecraft 
(P„y ) , when subtracted from Che landmark position in some coordinate 
frm%, will provide the measurement vector (M) 


-I “ ^I^E "^/C, 
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IcosAVsinAH 


U ^sinAM 


-Es/C - r + 1-90-6 
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the landmark being used » f (L,X«A^), 

the altitude of the landmark above the 

mean radius of the earth* 

longitude of the landmark* 

latitude of the landmark* 

vector position of the landmark relative 

to the center of the earth* 

vector position of the spacecraft 

relative to the center of the earth* 

altitude of the spacecraft above the 

mean radius of the earth* 

measurement vector from the spacecraft 

to the landmark* 

unit vector along M* 

the landmark tracker horizontal plane 

angular deflection from the boresight 

axis* 

the landmark tracker vertical plane 
angular deflection from the boresight 
axis. 


Figux^Q III-l 

Landmark Traaker Geometry - 
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Figure III-2 
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Accounting for hardware Disalignments, the same measurement vector in 
landmark tracker coordinate is 


h - 1^1 


" ^/C 


I 


“ l^E ^ “ 1^1 ^/Cj* 


From examination of Figure 111-2, the unit measurement vector in land- 
mark tracker coordinates is 
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sin AV 

cos AV sin AU 
cos AV cos AU 


However, the tracker instrument has no sensitivity to projections along 
its boresight axis* Therefore the tracker response to the unit vector 
U t will be 
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sin 
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UAy 

m 

cos 

AV sin AU 


producing a AH and AV as shown in Figure 111-2 as sensor outputs. 
Since Che sensor output will be corrupted by bias and noise, the sensed 
measurement will be 


AU ' 


s 


AV 

s 



where 


AH + bjj + Vjj 

AV + b„ + V 
V V 


^Hi component landmark tracker bias^ 

VyB component landmark tracker zero mean random 
noise, N (0, o 2), 

The component biases and standard deviations (o) are user selectables* 

The landmark tracker measurement may be compensated for knowledge of 
instrument bias. The bias knowledge may be a priori or through estima* 
tion. The compensated sensor output will be 
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Global Po»itionin£ SvitOTi ( CPS ) Model ^v«lopm«nt - The function of the 
received is to yield an eccinate of the uaar'a position end velo- 
city in inertial apace baaed on aignala from four GPS satellites* The 
detail associated with the nodeling the navigation signals from every 
GPS satellite is unnecessary for this siasilation. A aufficiant model 
simply adds measurement noise in the form of a position-velocity uncer- 
tainty to the user's true inertial position-velocity vector. This un- 
certainty is assumed to be made up of component uncertainties along the 
inertial coordinate axes. These components are assumed to be normally 
distributed with aero mean and standard deviation i ■ Xi, 

2l, Xi, Yi, Zi. 


The GPS measurement vector is therefore 

r X 1 





+ V + B(t) 


where 


Ic ■ 

X 

s 
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B(t) 


the CPS measurement vector (6 x 1), 

■ the user's inertial position vector (3 x 1), 

* Che user's ineritial velocity vector (3 x 1), 

{V 

* The position-velocity uncertainty vector X* 
components of which are ( (Q, o^^) ) (6 x 1) 

■ The position-velocity bias vector at time t. 



the 


When compensation is made for bias knowledge by subtracting that know- 
ledge, the resulting compensated measurement is 


- Z 

“8 -Z 


B 


where 

(B) * The position-velocity bias knowledge (b x 1). 

c. Star Tracker Measurement Model Development - When obtaining star 
Cracker measurements, it is important to identity when the fleid of 
view is blocked by a major body such as Che sun, moon, or earth. The 
following paragraph describes this process. 

Occultation - Occulation occurs when any major body comes between Che 
spacecraft and the target star or if a sufficiently bright object falls 
within the field of view of the tracker to render the device inoper- 
able. For Che simplified sC:*r tracker model, occultation will be de- 
fined to have occurred if the cracker bores ighc axis Couches any part 
of one of three major bodies— sun, earth, or moon. 
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Accounting for hardware misaligninanca, cha aana aMaauranent vacCor in 
landmark tracker coordinate is 




1^1 ” ^/Cj 


“ l^B ‘ 1^1 ^/Cj* 


From examination of Figure Ill~2, the unit measurement vector in land- 
mark tracker coordinates is 


Mi 

"wr 


Uixl Fsln AV 

Uiy • cos AV sin AH . 

Uiz cos AV cos AH 


However, the tracker instrument has no sensitivity to projections along 
its boresight axis. Therefore the tracker response to the unit vector 
U i will be 



UAx' 


'sin AV 

Ui » 


■ 



UAy_ 


cos AV sin AH 


producing a AH and AV as shown in Figure 111-2 as sensor outputs. 
Since the sensor output will be corrupted by bias and noise, the sensed 
measurement will be 



AH ' 

s 


-AH + bH^-V^ 

Zl - 


» 



AV 


AV + b„ + V 


8 


V v_ 


where 


W* component landmark tracker bias, 

Vy* component landmark tracker zero mean random 
noise, N (0, o2). 

The component biases and standard deviations (o) are user selectables. 

The landmark tracker measurement may be compensated for knowledge of 
instrument bias. The bias knowledge may be a priori or through estima- 
tion. The compensated sensor output will be 


AH 


AH - b„ 

c 


8 U 


s 

A 

AV 


AV - b„ 



La V J 
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Global Po»itionin£ Systtm (CPS) Mod«l D€Vtlopu»nt - Tha (unccion of the 
<«^S recaivad ia to yiald an astifflace of the uaar'a poaicion and velo- 
city in inertial apace baaed on aignala from four GPS aatallitaa. Tha 
detail aaaociatad with tha modeling tha navigation aignala from every 
CPS aatallita ia unnacasaary for thia aimulation. A aufficiant modal 
aimply adda meaaurenant noiaa in tha form of a poaition-valocity uncer- 
tainty to tha uaar'a true inertial poaition-valocity vector. Thia un-* 
certainty ia aaaumed to be made up of component uncartaintiaa along tha 
inertial coordinate axaa. Thaaa componanta are aaaumed to be normally 
diatributed with aero mean and atandard deviation i ■ xx, 

Zl, Xi, Yi, Zi. 


The CPS meaaurement vector ia therefore 




where 



+ V + B(t) 


^ ■ the GPS meaaurement vector (6 x Oi 

X • the uaer'a inertial poaition vector (i x 1), 

X - the uaer'a ineritial velocity vector (3 x 1), 

{V V*}^ 

y ■ The poaition-velocity uncertainty vector n* X * the 
~ componenta of which are ( (Q, o^^) ) (6 x 1), 

B(t) - The poaition-velocity biaa vector at time t. 

When compensation is made for bias knowledge by subtracting that know- 
ledge, the resulting compensated measurement ia 


Z • Z - B 
-g -g - 


where 

(B) ■ The position-velocity bias knowledge (b x 1). 

c. Star Tracker Meaaurement Model Development - When obtaining star 
tracker measurements, it is important to identity when the field of 
view is blocked by a major body such as the sun, moon, or earth. The 
following paragraph describes this process. 

Occultation - Occulation occurs when any major body comes between the 
spacecraft and the target star or if a sufficiently bright object falls 
within the field of view of the tracker to render the device inoper- 
able. For the simplified star cracker model, occultation will be de- 
fined to have occurred if the tracker bores ight axis touches any part 
of one of three major bodies— aun, earth, or moon. 
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Contidering Figure 

‘ ■ “"■* 1^’ '^1 ■ V<v*«>' ♦ <w>* * 


Since 

U, • &1 • • |U,I IR^I cos 0 . 


O ■ C08“^ 



where «nd Ug are known in the same coordinate syatem, and occul* 
tation ia denned to have occurred if a ^ 6. 

Under these circumstances, a star tracker update is not available from 
the occulted tracker. 

Konoccultation Star Tracker Itodel (Simplified) - If occultation has not 
occurred, a reference star will be defined to have been sighted on the 
star tracker boresight aria. In this simplified case, this will be 
modeled as a unit vector along the boresight axis. A bias and zero 
mean random noise will be added to the measurement vector. Therefore 
the corrupted unit vector in the i<^h g^gr tracker coordinates will be 





fb 1 


"v 


0 


X 


X 

M-U,+b.+v- 
— ai -al — al — ai 

0 


"y 

4- 

V 

y 


1 


_0 _ 


0 . 


The bias elements b,^ and by will be user selectable: The noise 

elements Vg and Vy will be zero mean normal variables N(0* a^) » 

where oVg and oVy will be user selectable. ^ 


II1-7 



Occulting Geometry 


Simulation Geometry 


Rp , Rni» Rg * line-of-sight vectors from the spacecraft (S/C) to the 
center of the earth, moon, and sun, respectively, 

s general radius vector to the 1^^ celestial object, 

* the radius of the 1^^ planet, 

* the unit vector to the target star 

« The unit vector along the star tracker boreslght axis 
(simplified model) , 

« {0,0,1}^ In tracker coordinates (simplified model) 

B « half-cone angle from the spacecraft to the 1^^ major 
body center and body limb, 

o ■ angle from the spacecraft to the 1^^ major body center 
and the target star. 

Figure III- 5 Star Tracker Meaeuremnt Model Development 
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Sinc« the star tracker has no sensitivity to projections along its 
boresighti only components perpendicular to the boresight axis in the 
tracker coordinates will be sensed. The sensed measurement for the 
i^h scar tracker will be 


AH 


bx + Vx 


m 


®i 


by + Vy 


The sensed measurement may be compensated for known biases. This bias 
is assumed to be known a priori and not estimated at this time. The 
compensated measurement vector is therefore 



&Vs^ - b^ 


Gyro Model - The general function of the gyro model ia to generate and 
process the angular changes sensed by the strapdown gyros at evenly 
spaced time intervals determined by one scan line of the science sensor* 

Processing consists of correction for gyro-considered parameters — 
drift, nonorthogonality, scale factor^ and misalignment* The compen- 
sated gyro output is used to update the quaternion so the attitude of 
the spacecraft can be known as accurately as possible at all times. It 
is then used to propagate the attitude state transition matrix, which 
is used in turn to propagate the covariance matrix* 

The major input to the gyro model is the angular changes of the space- 
craft during each sample period* The length of this period defaults to 
0*1 second. The input data consist of three variables that correspond 
to the angular changes about the three orthogonal body axes* 

Other input data required by the gyro model are the initial value of 
the attitude quaternion and the time at which it is valid (the last at- 
titude reference measurement time)* The attitude quaternion, which is 
propagated by the gyro model, is later updated by the Kalman filter 
when an attitude reference measurement is processed. 

The gyro model first simulates the gyro output for a specified inter- 
val. This is a three-step process in which the trif« attitude rate is 
obtained from an attitude profile table. These rates are transformed 
into individual gyro coordinates through 
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where 


• change in spacecraft attitude in body coordinates , 

>• transformation fr<Ms body to gyro coordinates accounting for 
misalignment of gyro package (the inverse of this transforma- 
tion is expanded below), 

■ transformation from gyro package coordinates to individual 
gyro coordinates accounting for nonorthogonality and scale 
factor, (the inverse of this transformation is expanded below), 

A6- ■ change in angle sensed by individual gyros. 

The change in angle sensed by the individual gyros is corrupted with 
gyro drift and random noise to yield the modeled gyro output 


A0„ - (A0^ + 8 + n) X AT 
M G 

where 

output, 

6 s gyro bias, 

T) » gyi'o noise 

AT ^ sampling interval. 

The gyro output is processed to eliminate gyro systematic errors (proc- 
ess noise). The systematic errors include the constant drift rate, the 
scale factor, the nonorthogonality of the gyro input axes, and the mis- 
alignment of the ideal gyro input axes with respect to the spacecraft 
axis system. These compensated angular chantes are calculated by the 
equation 


’ao; 


* “dx • 

A0^ 


A0, + ‘A*" 

0 

B G G G' 

2m Dy 



A0-, + • A^ 

3^ 


_ 3ro Dz J 


where 

bI'g * misalignment transformation matrix, 

GTq s gyro scale factor and nonorthogonality transition matrix, 
AO^, AQ^, AO^ > compensated gyro data (radians), 

AOijjj, ^® 3 in “ input data from gyros (radians). 
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^Dx» Woy, Wdj* gyro constant drift rate (rad/sec), 
At • gyro measurement interval (s). 


The matrix gTc has the following components: 


T > 
G G 



• 6 


1 


S 

y 


0 



s 

z 


where 

^xt Sy, $2'* gyi'O scale factors (radians /count) 


<Sl> <S2> 63 *° gyi**) misalignment angles (radians) 

Notice that G^G' is the inverse of the transformation used to model 6 
the individual gyro outputs. However, different values of S and 6 are 
used for modeling and compensation for realistic modeling of system un- 
certainties. 


The matrix B^G has the following components 


C4j • S4j 

“2 • “X 


where the abbreviations C and S are used for cosine and sine^ respec- 
tively, and A3, A2, Ai are Euler angles representing rotations 
about the z axis, y axis, and x axis, in that order* Again, B^G is 
the inverse of the transformation used to model gyro output but the 
values of A3, A2, and Ax are different. 

The compensated gyro data are then used to compute the differential of 
the quaternion that is used to propagate the quaternion. The quater- 
nion is used to update the attitude transition matrix* The net result 
of this process is an update in the attitude estimate and a new atti- 
tude transition matrix* The updated quaternion may be found through 


CA^ • 

CA, 




SA3 . 

CA2 




SA^ • 

SA2 • 

SA3 

-CAi 

• SA3 

SA3 • 

SA2 


+ CAj^ 

* CA 

CA^ . 

SA, • 

CA3 

+ SAj^ 

. SA3 

CAj^ • 

SA2 


-SA, 

• CA 
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®0. ei, 62, 63 

* old quaternion, 


eg. ej. e-. e' 

■ updated quaternion, 


A01, A02. A03 

* compensated gyro output. 


An overview of the measurement models is illustrated 

in Figure 1II-4 

%m No«e 
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Star TracAer 
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V 
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Sun — 
Moon — 


GPS Moaturomont Model 


Figure III-4 Ovevview of Mcaaux^.nent l-hdela 
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Dynamics Model and Navigation St<f Infgrator 

Dynamics Model - The dynamics model calculates the derivative of 
the spacecraft navigational state, which will be integrated to produce 
the navigational state vector. This is done in part by calculating the 
total acceleration of the spacecraft due to solar pressure and gravi- 
tation effects of the sun, moon, and earth, including fourth aonal har- 
monic terms. 

The sequential process for determing the forces acting on the space- 
craft is: 

1) Calculate Julian data 

2) Calculate solar and lunar positions 

3) Calculate solar pressure acceleration 

4) Calculate gravitational acceleration due to sun and moon 

3) Calculate gravitational acceleration due to earth using four aonal 

harmonic terms 

6) Calculate total acceleration 

Algorithm and Process - The position of the spacecraft is calculated by 
solving three simultaneous, second-order differential equations 


Xi “ -Xi‘ + 8l(t,X) + ai(t.X) 

X2 “ -X?* p + 82(t,X) + a2(t,X) 

X 3 “ “Xa* ^ + 83(t.X) + a 3 (t,X) 

where 

X -[Xi, X 2 . X3,]^ 

y “ earth gravitational constant (3.985491204E + 14 

R =• (Xt + x; + X§)^. 

Xj, X 2 , X 3 ■* coordinates of spacecraft, 

81* 82* 83 ” accelerations caused by zonal harmonics of earth gravity, 

^1* A 2 , a 3 » solar radiation pressure perturbations, sun and moon gravity. 


The equations for computation of the zonal harmonic accelerations are 

Xi 

81 " ^'l 

R3 


X 2 

g2 - -y-r Fi 

R^ 


83 -y — 

r2 


X3 

^ F, - 
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where 


4 


Fi - 


F2 - 




R 3 


R. ^ 

f2J2 

e 

R 


— • £ 4 ^ 4 * 


R ^ 



flJ2 + 

e 

R 

• f2>73 ^ 

r • ^ 3 ^ 4 . 


f 2 = -7.5 + 1.5. 

f 3 - -17.5 + 7.5 (^) , 

£4 = -39.375 (^)** + 26.25 - 1.875, 

J 2 *■ 1082. 7E -6 (Harmonic term in earth gravity model). 
J 3 » -2.56E -6 (Harmonic term in earth gravity model). 
J 4 » -1.58E -6 (Harmonic term in earth gravity model). 


b. Navigation State Integrator - The position state is advanced in 
time by numerical integration of the equations of motion. The second- 
order equations of motion are composed of the external forces acting on 
the spacecraft. The external forces consist of geopotential , lunar and 
solar gravitation, and radiation pressure. The integration routine is 
used to propagate both the actual state and the estimated state between 
measurement time. 

Algorithm and Process - A study showed that the Runge Kutta Gill (RKG) 
fourth order numerical integration method is optimal for this applica- 
tion. It is self-starting, handles variable step sizes, and is suf- 
ficiently accurate. The RKG method for numerically integrating dif- 
ferential equations is described here. 

The change in the value of the function during the computing interval 
is calculated by 

Ay - |(ki + 2 ( 1 -U)k 2 + 2 ( 1 +U)k 3 + k 4 ] 
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whera 


ki - Yn), M ■ »^/2 , 

k2 • h*f(t„ + ^sh, Yu + >jki), 

ka ■ h*f(tn + *sh, + y)ki + (1 - w)k 2 ), 

k4 - h*f(t^ ■*“ **» yn ’ 

h * computing interval (seconds)* 

tQ ■ time of beginning of computing interval (seconds) , 

> value of function at beginning of computing Interval. 


The derivative function f shall be evaluated four times to calculate 
the change in the function being integrated during the computing inter- 
val. 

At first appearance, the software algorithm shown in Figure 111-5 has 
no relation to the mathematical description of the integrator. The 
following is intended to show that the two are indeed identical. 



Initial Entry 


Null 0 Matrix 


Compute Position Dynamics 



• 1 or 3 


' - T + DT/2 


TP - AA(K) X [ZD(D - BB(K) x 0(1)] 


W - Z(I) 


Z(I) - Z(I) + TP X DT 


TP - [Z(D - W]/DT 


0(1) - 0(1) + 3 X TP - CC(K) X ZD(I) 


where AA - 0.5, 1 - 2/2, 1 + 2/2, 0.5/3; BB 

2/2, 1 + 2/2, 0.5 


Do Until I ■ 6 


Do Until K - 4 


Figure III-5 SoftMore Algorithm 




- 2, 1,1, 2; CC - 0.5, 1 - 
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Proceeding through Che outer "do" loop for K ■ 1 Ce 4, we find the fol 
lowing results 


k " 1 




ti - Hf 

yi ■ yo + HKitf 

qi - f 


k - 2 


f 


t2 

y2 


(1 - (f - f) 




yo + % + dt(l - (f • - f) 


yo + dt(-5j + ^) f + dtd - ^)f 


<12 


f + 3((1 - -y) (f - f)) - (1 - y)f • 


Jl /2 

f - 3(1 - y) f + 2(1 - y)f ’ 


f" 


f(tn + y. y + dt(-»s + y)f + dtd - y)f 


t 3 - (1 + y) (f" - f + 3(1 - y)f - 2(1 - y)f ') 


ya • yo + dt(-Js + y)f + dtd - y)f + dt(l + y)(f" - f + 3(1 - 
*^)f - 2(1 - y)f •) • yo - dtyf + dtd + y)f" 


2 /* 2^* ^ . 2 
qj - f - 3(1 - y)f + 2(1 - y)f + 3((1 + y){f" - f + 3(1 - y) 


2(1 - - (1 +y)f") 


2'‘ ^ • 2 
~ (1 2 ) ^ 2(1 + 

f” - f(t„ + dt. yo - dt&' + dtd + •y)f") 


ti» - l/6(f" - 2 - Ijf - (1 + y)f + 2(1 + y)f") 


l/6(f" + f + 2(1 + y)f - 4(1 + y)f") 


y4 “ yo - dtyf + dtd + y)f" + l/6(dtf + 2dt(l + y)f - 
4dt(l + y)f" + dtf ”) 

- yo + l/6(dtf + 2dt(l - y)f + 2dt(l + y)f" + dtf'")* 


lll-lb 



Making Che subaCiCuCiona 
kj ■ dtf 
k2 - dtf 
k$ ■ dtf" 
k4 - dtf" 



yielda 

y - yo ■ Ay ■ l/6{ki + 2(1 - k)k 2 + 2(1 + k)k 3 + k 4 ), 

which ia identical to the original algorithm. 

The Extended Carlaon Square Root Filter Mechanization 

The Carlson square root filter is an extension of Che conventional Kal*- 
man filter where the optimal gain, state and covariance updates are de- 
termined using a matrix that is the square root of Che conventional co- 
variance matrix. The advantages of this mechanization are higher ac- 
curacy for a given machine precision and a guaranteed positive covar- 
iance. The extended filter is necessary due to Che nonlinear nature of 
Che process. The extended filter propagates Che state, state transi- 
tion, and process noise covariance by integration of differential equa- 
tions rather than linear propagation over the appropriate intervals. 

The model equations used for the filter mechanization are the state 
differencial equations 

X « f(X (t), t) + W(t) , W(t) * N(0, s(t)) 

where 

Q(t) - E(W(t)W(t)T), 

E ■ expectation operator, 

W(t) B normal random variable process noise with zero mean and co- 
variance, Q(t), 

f * a nonlinear function of X(c) and C, 
the measurement equation 

ik ° hk(X(tk)) + Vji , Vk • N(0, Rk) 
where 
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R|i * the covariance of Vk 


Vk ■ norinal random variable measurement noiae with aero mean and co- 
variance, Rk, 

and the measurement estimate equation 

ik • llk(X Ctk) ) . 

The following definitions are made 


■ i 

H(t) - H 


X - X(ti) 
X - X(ti>. 


The following nonlinear differential equations are integrated by 
Runge-Kutta Gill fourth order numerical integration over the interval 
(to, ti). 

State Estimate 

X - f (X(t) , t) X(to) - Xo 

where f(x(t), t) Is the dynamic model 


State Transition Matrix 

- F(t)*(t, to). ♦(to. to) - I 

Process Noise Covariance 

- F(t)Q(t) + Q(t)F(t) + Qo . Q(to) - 0. 

Once the state transition matrix differential equation is integrated 
over the interval (tQ, t^), the state covariance may be propagated 
forward to t^ by the equation 


The measurement residual is defined as 
P(ti) - ♦(ti, to)P(to)^^(ti. to) +Q(ti). 

where 

^(tj) » the respective sensor measurement model output at time ti 
Z(ti) « Z(ti) - i(ti) - Z(ti) - h(X(ti)) 
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Under normal Kalman filter mechaniaations, the optimal gain, atate aa- 
timate update, and atate eatimate covariance update equatioaa at the 
measurement time are reapectively 

k(ti) - P(ti)H(ti)(H(ti)P(t,)H'T(ti) R(ti))“^ 

X*'(ti) - X(ti) + k(ti)i(tj) 

P'*‘(ti) - (I - k(ti)H(ti))P(ti) 

However, with the Carlaon aquare root mechanisation the aquare root of 
the state estimate covariance P(tj^) is found by Cholesky decomposi** 
tion so 

S(t^)S(t^)^ - P(t^) 

SCt]^) is then used along with H(ti) to find the optimal gain and 
state estimate update. The state estimated covariance update is found 
in the square root form S+itj) and the reconstruction is 
implemented by 

This process required the restrictions that measurement vectors be de- 
composed to scalar form and the measurement noise covariance be diagon- 
al (uncorrelated elements). The latter appears to be true for our mea- 
surement schemes and coordinate fraiaea used. The former is mechanised 
by using the elements of the measurement vector Z(t^) in sequence 
without advancing time from t^. This requires the use of the ele- 
ments of H(t|) by rows. 

The primary function of the remainder of this chapter is to develop the 
expressions for H(t^) for each measurement sensor that may be used in 
the filter. 

a. General Formulation 


The estimation state vector is 

y* *t y» *• ^0* ^1* ^2* ^ 3 ^ • 

The measurement model for each instrument has been previously described 
and has the form 

Zj • hj(X) + Vj, J - G, S. L 
where 

G ■ GPS measurement, 

S ■ Star tracker measurement, 

L ■ Landmark tracker measurement. 

The estimated measurement will therefore have the form 

Zj - hj(X) 
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where 


Z end X ere escimeces of the respective vector veriebles. 

Therefore 

dZj 

H4 “ “T" 

3X 

where 

Zj es e function of X hes been developed ee previously discussed, 
b. GPS Meesurenent Distribution Metrix - 




LITH 


h ■ hoQt) 


I 0 I 0- . ! 

X) 1 — .j X • * . 

0 I I I 0' « 

L I • J V 


Therefore. 

32 

H. • — - 


32 fl I 0 I O'] 

, » — ■ I—}— - 4 — --It where 

' ax |o } I I O'J 

[ i 0 o“| To 0 o"| To 0 0 dl 

0 1 0 . 0-0 0 0 . and 0 ' - I 0 0 0 o| . 

001 000 [_ooooj 


Star Tracker Measurement Distribution Matrix - 




Z. - h,® + y 


Z - h (X) 
—8 8 “ 


AH 


8 

m 

AV 


8 



Fin*HOs^ ^/cos 
sin” 


8(8in“^Us^))l f Ui 
-«»x J'l°- 


tihere Usx and Usy are elements of an estimated unit vector to e sighted 
star. Us. in the coordinates of the tracker of interest* (For small 
angles AHg « Osy and AV'a « Usx) 
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To sinulate cho measurement process ss closely ss possible, s unit vec> 
tor to a star is first established in inertial space# This would be 
the normal star measurement. This unit vector is established using ex* 
act star tracker alignment and exact body orientation and is then 
transformed into tracker coordinates using estimated parameters to de- 
velop the estimated measurement 


iia 


m 


0 

0 

1 


Ideal unit vector along the S/T boresight in S/T coordinates 


Transformation to inertial apace 

- iTs Us • iTg bTs Os 
where 

ITb ■ exact body to inertial transformation 

* ^iCqOi qit q2t < 13 ) (<li * body quaternian elements) 
bTs ■ exact S/T to body transformation 

• ^2^®» Y,{0}) (o, $ , Y ■ misalignment elements 

(0) ■ nominal orientation array). 

U| is the simulated unit vector to the star using exact attitude in- 
*formation. Now transforming this vector to star tracker coordinates 
using the estimated attitude information yields: 

Us ■ sTi Ul • sTb BTi Ui - sTb B^I iTb bTs Us 

where 

BTi ■ f3(q0. qi. q2. q3). (qi ■ estimates of qi), 

sTb ■ f4(a, 6i Y, {0}), (a, 6,Y •• estimates of o, S, Y ) 

where 

a, 6, and Y are a priori estimates of misalignment angles of the star 
tracker of interest relative to the nominal alignment. These parame- 
ters may be estimated at a later data but will be used here as implicit 
consider parameters. 

{0} is a nominal orientation array. 
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Since Zg is a nonlinear function of the elementa of the state eati** 
mate X, Ha otust be found by 



ddHs* 

3 X 

3 AVs 

A 

ax 


s 


'ausy 

a^ 

aOsx 

ax 


i 8ln“^U8X 

ax 


Uhere 

aAVs ^ 1 ausx s au8x 

" A-(us^^ avi 


a^ 

^xi 




/ 6 »y _ Y 

\(:os(sin"^08x)/ 


aUsy ^ ^ 1 aUsx 

C 08 ( 8 ln~^Usx)axi + Usy^Us x */ l-fesx)^ 3 x 1 
(cos(sin“^Usx))^ 



** y» 2^ Z} <i0^ , q2» q3 • 


Since ^73 is independent of any elements of the estimated state, 
consider 


l'*^s “ l'*^B b'*^S 


®n 

*12 

*13 

821 

*22 

*23 

*31 

*32 

*33 
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Therefore 


®13 <lll <I12 fll3 

U8 ■ gfg jfj 823 " S^B ^21 922 923 

833 931 932 933 


where 


9n ■ 9o 9r “ 92 - 93^ 

912 “ 2(9192 + 939o) 

913 ■ 2(9193 - 929o) 

921 “ 2(9192 - 939o) 

922 “ 9o^ “ 9l^ + 92^ ” 93^ 

923 ■ 2(q293 9l90> 

931 ■ 2(9193 + 929o) 

932 ■ 2(q293 “ 9l9o) 

A ao ^2 ^ 2 ^ ^ 2 

933 ■ 90 “ 9l - 92 + 93 


Since a I |S , end V are not being estimated in the filter, they are con 
sidered constant and imbedded in gig as 

bii bi2 bi3 

8^8“ **21 b22 b23 ■ * 

b31 b32 b33 
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There fore 



ai3 


Osx 

a23 

• 

dsy 

«33 


Osz 

«» m 


Since Usz is not observed by the star tracker, it is dropped from the 
measurement vector* 

Usx ■ bn (qnai3 + qi2«23 ^13«33) 

+ bi2 (^21*13 ^22*23+ ^23*33) 

+ bi3(q31*l3 + q32023 + ^33*33) 

flsy • b2i (qiiai3 + qx2«23 + qi3a33) 

+ b22 (q2iai3 + q22«23 + q23«33) 

+ b23(q3iai3 + q32«23 + q33®33) 


Define 


M - [fe] 


ausx 

3x1 

dOsx 

Iqi 


ausy 

3qi 


•^1?^ « 0 for xi » X, y, z, x, y, i 
3x1 

^u(aX 3 l--a 3 %- 33 %) ^ 

+ bi 2 (ai 3 T^ + a 23 ^ + a 33 ^) 

+ bl 3 (ai 3 ^ + * 23 ^ ® 33 ^) 

‘ b 2 l(ai 3 ^^ + ® 23 ^^ ■*“ * 33 t|^) 

+ b 22 (ai 3 %: + a 23 ||f- + * 33 ^) 

+ b 23 (ai 3 ^fji + ® 23 ^|^ *34 
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where 


la«. . 2q. Ml* . 

2qa 

3a, a 
34o 

Mw. —2a, Ma* . 

3<io ^ 3^0 

24o 

3a., 

34o 

. 2qa - 

Ho 3qo 

2qi 

Mil 

3qo 

Or in 

alternate form 



3 f 

‘qij^ 


A 

‘lo 

^a 

A 

-‘12 

3^0 ‘ 

- 2 

-‘la 

A 

‘lo 

A 

*li 




.^2 

-^1 

A 

\ 

Similarly 





3 

3^1 

[qij 

■ ^ 

A 

*11 

^2 

A 

^2. 

A 

-^1 

A 

“la 

A 

‘Iq 




Us 

“^0 

■’l 

3 

3^2 

N 

- 2 

A 

-‘»2 

A 

A 

m 

A 

-’o 

A 

"a 




. ^0 

^a 

-«a 

m 

d 1 

■qlj] 


Pa 

-•la 

A 

‘lo 

A 

sqg 1 

« 2 

A 

-qo 

A 

-’3 

A 

’al 





’a 

‘a. 


•2^2 


2^1 


2qo 
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So 



Finally 


Hs 


’jAH 

9AH 

9AH 

3AH 

3AH 

3AH 

3AH 

3AH 

3AH 

3AH 

IT 


9z 

3x 

3y 

IT 


3qi 

3q2 

3q3 

dAV 

3AV 

3AV 

3AV 

3AV 

3AV 

3AV 

3AV 

3AV 

3AV 

- 

3y 

3z 

3x 

3y 

3 z 

3qo 

3qi 

3q2 

3qa 

1 

1 

1 

3AH 

3AH 

3AH 

3AH 





1 

1 

1 * 

3*10 

3qi 

3q2 

3qa 

. 




1 \ 
1 

1 

3AV 

3AV 

3AV 

3AV 





» 

« 

1 

« 

3qo 

3qi 

3q2 

3qa 






The approximation equation will be used initially and expanded later* 
d. Landmark Tracker Distribution Matrix - Hil 


Zl = hJl(X) + NH 
and 



AHL 

bJUX) 

a 


1 

\Vl 

m «M 

c: > 

1 





A 

sin-Huty/c 
_ sin~^ 


i*‘Huty/cos(sin~^U)lx)) 

A 

•^ULx 


A A 

where Utx and I'ty ^re elements of an estimated unit vector to the 
sighted landmark, Ut . in landmark tracker coordinates* (Again for 
small angles AHJl ^ Utx and AVt = ULy*) 


The simulated measurement is made by unitiaing the vector from the 
spacecraft to the Landmark. The result is then transformed to tracker 
coordinates through the estimate of the transformation from inertial 
space to landmark track coordinates, j^T^. 

The estimated measurement vector in inertial space if formed by 

Bl ■ (iTgiE) -il 

where 

W. « the landmark position vector in earth fixed coordinates , 

* 

Pj = the spacecraft position vector estimate in inertial space. 

Transforming this to the best knowledge of landmark tracker coordiantes 
yields : 


. gTi ((iTg I^) - li) 

* £^8 B^I ( ^I'^E 

where 

bTi * f3 (qo, qi, q 2 . q3^ 

Tfi = f5 (a, Y, {0}), (o, 0, Y a estimates of the tracker 

misalignment angles. tO} > 
nominal orientation array). 

Consider that 

iTe Le * Li “ l 

il « £ - fe (x, y, 1 ) 

gTi * [qij] (defined in star tracker section) 
jlTg » laij]. 

So 

UijJ IqijJ (L - P). 

A 

Unitizing M^into a unit vector 


U£ 



”u£x 

U£y 

_U£2 


T A L 

where |M£| - (M£ M£)^ . 


since Z is very similar to Zs, H may be found directly as 


Hi 


3^ 

3x 

‘3 

3AH1 

3x 

3AVt 

3x 


3x sin 

sin' 

3x 

cos(sln“^Utx)) 

"^Ulx 


“aOjty 

3x 

autx 

3x 


where 

3Ayt ^ I aptx aOtx 

3xi /l -(Uix)‘ ~ 


and 

3AHt 

3x1 


^ . -5- * . . _J. 

cos(sln ^Utx)3xj Uty'»Uty«Utx»/l - 

j\ / y.- ■ ^ (cos(sin~lUZy))^ 

y ycos(sin'lO£x)j J 

3Uty 

3x£ 


Xi = X, y, z, X, y, z, qQ, qj, q2, qa 


*i “ X, y, z, X, y, z, qg, qi. q2. q3* 
Since 


Vl 


ML. 

|M£| 

im| 


am Mt 9lE«| 
3>ci ~ “T^i 



. 3Utx 

iUislLiJJi 


3x1 


: a three vector for each xi 

making a total of (3x22) 


where 


am 

3jti 


{aij} {qij} (L-i) + {qiJ} ^ 


1 - 1 . 10 


^ - -il mJ)& + . (^T ^ ^ 


. (4- 

ax^ 3x1 


{qiJ} “ iO} for XI - X, y, z, i, f, I 

3x1 


*lo ^3 “^2 

"^3 ^0 ^1 

AAA 

q£ ~^i ^0 


for Xj - Mq 


A 

*J2 

5a 

A 

■q2 

A 

A 


<l2 

-qi 

A 

-<io 

A 

1l 

A 

«l2 

A 

*l3 


+q3 

A 

% 

A 

-qi 

A 

"<Jo 

^3 

A 

-<12 


for • qj 


for XI • q^ 


» 2 


and 


-q3 qo <11 

-qo “q3 ^2 

qi q2 qa 


for xi “ qj 


^ - i] = 0 for ii - K I. q^. qj, ^2* ^3 
for XI “ X 

A A 

for xj^ ■ y 

A A 

for xi ■ z # 


-1 

0 

0 

m I 

0 

-1 

0 


0 

0 

-1 
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B. 


EXTRACTED GCP MEASUREMENTS 


The second mode of operation allowed by GCPSlM provides for a mixture 
of modeled and actual measurements. The CPS and star tracker measure- 
ments are modeled as discussed previously. The GCP measurements , how- 
ever» are extracted from actual imagery obtained from Landsat. For 
this process 9 two passes of corrected Landsat imagery over the same 
area are required. The first pass is used to extract CCPs to be incor- 
porated into the GCP data base discussed in Chapter II. Since the ex- 
tracted GCPs are used for references^ it is desirable that the first 
pass of imagery be of good quality with as little cloud coverage as 
possible. It is also necessary that the latitude and longitude of the 
four corners of the image be known so the latitude and longitude of the 
individual landmarks can be calculated. Given these constraints, it is 
possible to construct a GCP data base to be used in the simulation. A 
program, entitled LMEX, has been developed to aid this process. 

The second mode of operation functions much the same way as the first 
mode. The sequencer defines the type and time of the star tracker and 
CPS measurements. The GCP data base containing the number of scan 
lines to the next GCP is C 0 ;isulted to see if the GCP will be encoun- 
tered prior to the measurement; if there is not enough time, the mea- 
surement defaults to a GCP sighting. The primary difference between 
the two modes of operation is that prior to exiting from the sequencer, 
an interrupt is sent via a parall«:l data bus to a separate task hosted 
on the PDP 11/70. The interrupt contains the number of lines to be 
scanned prior to measurement time and a flag identifying whether a GCP 
sighting will be available or not. 

The entire purpose of the second program is to perform all the func- 
tions of the registration processor and to simulate the sensor front 
end. 

The program has been separated into three tasks running independently 
on the PDP 11/70 at different priorities. Intertask communication is 
performed through a shareable global area, and synchronization is con- 
trolled using global event flags. The three tasks are described as: 

1) SCRLL - Accepts an interrupt from the VAX 11/780 with the number 
of lines to be scrolled and a flag indicating whether a GCP 
sighting is available* SCRLL also performs the actual scrolling 
of imagery across a video monitor to simulate the acquisition of 
imagery; 

2) SCOUNT - When enabled by SCRLL between each scan line, SCOUNT up- 
dates the line counter, which indicates the distance to the 
search area. If the line counter is 0 or negative, SCOUNT also 
stores the search area imagery as the data are scanned; 


i 
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CORR - SCOUNT initializes CORR when the entire search area has 
be read into memory. CORR then perforins the SSDA registration 
process using the reference GCP extracted from the first pass of 
imagery. CORR also passes the registration vector back to SCRLLi 
which passes it to GCPSIN. 

The interrelationships of the tasks are illustrated in Figure Ill-'b and 
the process is illustrated in Figures For a detailed descrip- 

tion of the software) consuit Volume III* 

Several problems were encountered when developing this simulation. 
First) both passes of imagery must be corrected tor distortions because 
the distortions associated with the mult ispectral scanner (HSS) are too 
dynamic to be modeled. This problem will not be encountered with 
either the thematic mapper or the MLA because distortions are determin- 
istic. If uncorrected imagery is used) the location of the extracted 
CCPs cannot be accurately determined* and the registration vector will 
have a large component caused by this uncertainty and the distortion of 
the second pass. As a result* the filter has poor convergence charac- 
teristics because the effective noise* rather than being 0.1 pixel* is 
now on the order of 10 pixels. The imagery acquired as GFE was uncor- 
rected imagery and therefore could not be used. 

The second problem encountered is that the true vehicle state at the 
time of scene exposure is required. However* if the imagery has been 
corrected* this state is no longer valid because the errors associated 
with vehicle state have been extracted. 

In an operational system* these problems will not be encountered be- 
cause the distortions associated with the MLA are totally deterministic 
so the location of the GCPs are well-known functions of spacecraft 
state* 
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H)c GCP tiBulatlon tystMi configuration diaeuaalon in ebaptar II will 
not ba rapaatad bare. Otbar araaa aueb at aaauaad orbital ebaractar* 
iatieci vabicla notion, aanaor arror budgat and tba raaulta of para- 
natrie atudiea will ba diacuaaad in detail* 

For tba purpoaa of tbia atudj, tba apacacraft orbit waa aaaunad to ba 
circular with an inclination to tba aguatorial plana at 80 dagraaa* 

For aaaa of aatup and intarpratation of raaulta, tba beginning of aacb 
ainulation run waa taken to ba at a poaitiwa aacanding node wi^ tba 
inertial coordinate fraan coincident with tba aartb'a coordinate 
fraoM. Tbia did not detract fron tba generality oi tba raaulta ainca 
apacific acenaa on tba aurfaca or apacific atar cataloga ware not uaad 
in tba ainulation. 

Tba vabielc'a notion waa initially intended to ba tbat of table lookup 
of vabicla body ratca about aacb vabicla axia, Jollowad by integration 
of thaaa ratca to datamina attituda* Tbia, bowavar, waa ebangad to a 
forcaa nadir-oriantad vabicla through derivation of a vabicla body rata 
vector found by conputing tba croaa product of tba true vabicla veloci- 
ty and poaitlon vactora* Tbia rata vector waa uaad aa input to both 
tba gyro nodal and tba trua attituda propagator wbara true attitude waa 
obtained through integration of tba body rates* 


ERROR BUDGET 


I 


f 

I 




1. Attituda Rafaranca Unit 

Tba attituda rafaranca unit for tbia ainulation waa aaaunad to ba of 
tba NASA atandard type (DRIDU II)* It conaiata of three noninally or- 
tbogonal gyroa wboae cbaractariatica are fairly wall docunantad. lv<a 
apacificationa for tbaaa gyroa ware atudiad and the tenn applicable to 
tba aaaunad ayatan configuration and anviroonant ware uaad* The para- 
natara apacifiad and uaad are outlined in Table IV-1* 

2. Star Trackara 


Each of the two atar trackara aaaunad to ba on board tba vabicla are of 
tba NASA atandard type. For tba purpoaaa of tbia ainulation, tba opti- 
cal axaa of tba two davicaa ware aaauawd to ba orthogonally oriented in 
a plana containing tba vabicla poaition vector and perpendicular to tba 
vabiclt' velocity vector. Docunantation, including tba 0AD8 and axiat- 
ing apacificationa, waa atudiad and tba arror tama fait to ba applic- 
able to tbia ayatan configuration and environment ware uaad* Tbaaa 
paraMtera are outlined in Table IV-2* 
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Table IV~1 Attitude Referenoe Unit Error Coefficients 


ParMMter 

Value of 
Conven- 
tional 
Units, lo 

Value Used, lo 

Source 

Gyro Long-Term Bias 

0.017 deg/b 

0*08 X 10*8 tad/t 

Spec 

Uneertalnity 

Scale Factor Uncertainty 

90 ppm 

90.0 X 10“0 

Spec 

Monortbogonality of Mounting 

2 arc-a 

10 X 10"^ rad 

Assumed 

Package Alignment Uncertainty 
Selative to Mounting Surface 
Bach Axis 

10 arc-a 

50 X 10-6 

Assumed 

Gyro Noise Equivalent Angle 

1*3 arc-s 

10 X 10-6 rad/s 

Spec 

Absolute Bias 

2 deg/h 

10 X 10-6 rad/s 

Spec 

Absolute Scale Factor 

0.5X 

0.005 

Spec 


Table IV~2 Star Tracker Error Coefficients 


Parameter 

Value of 
Conventional 
Units, 1 0 

Value Used, lo 

Source 

Bias, Bach Axis 

10 arc-s 

50 X 10-6 rad 

Spec 

Quantisation 

1 arc-s 

5 X 10-6 rad 

OADS 

Misalignment 
Bach Axis 

5 arc-s 

25 X 10-6 rad 

OADS 

Star Catalog 

1 arc-s 

5 X 10-6 rad 

OADS 


3. Global Positioning System (CPS) 

The GPS is assumed to be available to the spacecraft virtually concinu 
ously. The update interval using GPS was predominantly i seconds; how 
ever, variations have been made in this paraiseter to study its affect 
on system performance. The GPS system error parameter assumed for the 
purposes of this simulation are presented in Table lV-3. 


Table IV~3 GPS Error Coefficients 


Parameter 

Value of 
Conventional 
Units, lo 

Value Used, 1 o 

Source 

Velocity Bias 
Velocity Moisr 
Position Bias 
Position Noise 

0 

0.005 m/s 
0.04 m 
4.0 m 

0 

6 X 10-6 |ai/« 
4 X 10-5 las 
4 X 10-3 )qi 

Assumed 

Assumed 

Assumed 

Assumed 
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Landmark Trackar (LMT) 


Tba UfT was aaauaad co have its optical axis collnaar with tbe nadir 
pointing axia of tbe apacecraft. Tbe look angle to tbe landmark ia in 
tbe plane containing tbe optical axist perpendicular co tbe vebicle 
velocity vector. It baa a maximum value of 10 degreea with two modea 
of operation. Tbeae aK>dea are uaer-aelec table to take on either fixed 
a^cified look anglea or unifocnly diatributed look anglea of a apeci- 
fied maximum magnitude. Table 1V>4 ouClinea tbe error cbaracteriatica 
aasumed for tbe LMT. 

Table IV-4 Landmark Tracker Error Coefficients 


Parasieter 

Value 

Conventional 
Units 1 a 

Value Used 1 o 

Source 

Bias, Each Axis 

2 arc-s 

10 X 10“6 rad 

Assumed 

Quantization 
Each Axis 

1 arc-s 

5 X 10“6 rad 

Assumed 

Misalignment 
(three axes) 

5 arc-s 

25 X 10~6 rad 

Assumed 

Landmark Location 
Uncertainty 

1 arc-s 

5 X 10”b rad 

Assuaied 


B. PARAMETRIC STUDIES 


The intent of tbe parametric studies was to determine the system's 
sensitivity to such parameters as LMT update interval, look angle mode 
and range, GPS update interval and accuracy, and recovery time after a 
long period without LMT updates with and without star tracker backup. 
The study was progressive in tbe above order using as parameters the 
values that appeared to be optimum from tbe previous study. Witb tbe 
exception of tbe study of system sensitivity to tbe GPS accuracy, the 
error budget referred to was used uniformly for all studies. 

Table IV-5 outlines the cases studied by title and aiost significant 
parameters. Plots showing tbe results of these runs have been included 
at tbe end of this chapter. 

1. Sensitivity to GCP Update Interval 

The sensitivity to GCP update interval was investigated by examining 
the system performance parameters of estimated position, velocity, and 
attitude accuracy after system convergence for GCP interval variation 
from S to 320 seconds in factors of 2. These are respectively cases 1 
through 7. While staking these variations, tbe GPS update interval was 
held at 5 seconds and tbe error budget held to tbe values previously 
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specified. Both random (♦lO-deg maximum) and fixed (0-deg) look angles 
were considered for the landmark tracker pointing for each case pre- 
sented. No star tracker measurements were used during this study. 


Table IV-6 Cases Studied 



GCP Update 

GPS Update 


Title 

Interval , s 

Interval, s 

Remarks 



GCP Update Interval Sensitivity 

CASE 1 

5 

5 

CASE 1 through CASE 7 were run 

CASE 2 

10 

5 

with -t-lO-deg random look angles 

CASE 3 

20 

5 

as weTl as a 0-degree fixed look 

CASE 4 

40 

5 

angle. 

CASE 5 

80 

5 


CASE 6 

160 

5 


CASE 7 

320 

5 


Look Angle Sensitivity 

C.SE 9 

80 

5 

0-deg fixed look angle. 

CASE 10 

80 

5 

2.5-deg fixed look angle. 

CASE 12 

80 

5 

3-deg fixed look angle. 

CASE 14 

80 

5 

10-deg fixed look angle. 

CASE 17 

80 

5 

■•-2. 3-deg random look angle. 

CASE 18 

80 

5 

■^3*deg random look angle. 

GPS Sensitivity Update Interval 

CASE 21 

80 

10 

CASE 21 through CASE 26 were 

CASE 22 

80 


run with ^10-deg look angles. 

CASE 23 

80 

40 


CASE 24 

80 



80 


Accuracy 

CASE 23 

80 

5 

GPS velocity error increased by 4X. 

CASE 26 

80 

5 

GPS position error increased by 4X. 

Recovery 

Sensitivity 



CASE 30 

80 

5 

^10-deg random look angles; linear 
cloud table; default to star 




trackers* 

CASE 31 

80 

5 

♦lO-deg random look angles; GCP 
shut down for 1000 s with default 




to star trackers. 

CASE 32 

80 

5 

♦10-deg random look angle; GCP 
shut down for 1000 s without 




default to star trackers. 
























The standard deviation of error is found by Identifying the peak at> 
titude and position errors after filter convergence* These values are 
taken as the 3 a attitude error » inplying that the error standard devia** 
tion a is the peak error valve divided by 3. Table IV**6 lists these 
standard deviation values. 


Table IV-6 Results of GCP Update Interval Sensitivity 


GCP Update 
Interval, s 

Peak Position 
Error n 

Peak Attitude 
Error, arc-s 

Bounds of Implied 
Surface Error, m lcr| 

5 

7.0 

16 


10 

6.5 

14 


20 

7.0 

11 


40 

7.0 

10 

12-14 

80 

6.5 

16 

18-20 

160 

9.0 

15 

17-20 

320 

6.0 

15 

17-19 


Initial runs showed that idien the GCP look angle was held at nadir* at- 
titude error convergence was poor. After some analysis, this phenOM- 
non was shown to be an observability problem. Uben uniformly distri- 
buted random look angles with peak magnitudes of 10 degrees were iaple- 
aiented, convergence became more rapid and definitive. For this reason 
the results were taken on the basis of random look angles. The results 
showing the attitude error convergence with look angles fixed at sero 
are shown for information only. The fixed look angle cases are indi- 
cated on the curves. 

In the way of quantitatively evaluating the results, the most desirable 
parameter to consider is the surface location error fiP resulting from 
the combined navigational position error 6P^ of the spacecraft and 
the implied surface position error 6 P 2 resulting from the attitude 
error je* former is printed directly by the computer program. The 

latter may be computed by 

6 P 2 = 705 km sin 66 ; 705 km 60. 

It may be shown that the position error on the surface 6P is bounded by 

{6Pi^ + 6P2^}^ 1 “SP ^ 6P1 + 6P2 

Table IV-6 shows the upper bound of the boresight determination error 
for each update interval. The peak navigational position error is 
fairly uniform at approximately 7 meters. The peak attitude error 
seems to have a minimum around the 40-second update interval point. 

This suggests a minimum position error implied at the surface to result 
from a GCP update interval near 40 seconds. The higher attitude errors 
for longer intervals seem to arise from the effects of gyro drift. 
Speculation on the high attitude errors for very short update intervals 
suggests the cause may be interference or rectification resulting from 
the closeness of the update intervals of GCP and GPS. Because the 
total state is not observable using landsurk tracker measurements 
alone, a short update interval causes artificially high covariance 
values. 
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Sensitivity to GCP Look Angles 
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Cases 9, 10, 12, and 14 examine the sensitivity of system performance 
to constant-amplitude look angle. The cases mentioned examine Che ef- 
fects of constant look angles of respectively 0, 2.5, 5, and 10 degrees 
with ocher conditions as previously specified and Che GCP update in- 
terval set at 80 seconds. Cases 17 and 18 examine Che sensitivity to 
look angles of uniform distribution between peaks of ^2.5 and +5.0 de- 
grees respectively again with an 80-second GCP update interval? All 
other conditions remained the same, and no star tracker measurements 
were used. 

The cases 9, 10, 12 and 14 showed an insensitivity to the magnitude of 
constant look angle with what appears to be a very long convergence 
time. As previously discussed, this is probably due to lack of total 
observability in this measurement situation. With random look angles 
convergence is more rapid with a +S-degree variation Chan with a *2.5- 
degree variation. However a comparison of case 12 with case 3 shows 
virtually identical convergence characteristics. The explanation is 
not apparent at this time. 

3. GPS Update Interval Sensitivity 

Cases 21, 22, 23, and 24 examine the system performance sensitivity to 
the GPS update interval. The intervals used were respectively 10, 20, 
40 and 80 seconds. The GCP look angle was maintained random at +10 de- 
grees maximum, and the update interval was set Co 80 seconds while all 
ocher parameters were held at their previous values. 

The peak values of error for position, velocity and attitude are pre- 
sented in Table IV-7. This Cable shows an apparent peak in position 
and attitude errors at or near the GPS update interval of 40 seconds, 
falling off on both sides. No explanation is readily available for 
this result except that possibly no 3a numbers were encountered due to 
the small number of data points. 


Table TV-7 Results of GPS Update Interval Sensitivity 


Case 

GPS Update 
Internal, s 

Peak 

Position 
Error, m 

Peak 

Attitude 
Error , s 

Implied Surface 
Error Bounds, m 

5 

5 

6.5 

16 

18-20 

21 

10 

7.0 

13 

15-17 

22 

20 

8.2 

16 

18-21 

23 

40 

12.0 

15 

18-21 

24 

80 

7.5 

13 

15-17 


1 





GPS Accuracy Sensitivit 


To detemioa bow sensitive tbe total system is to tbe accuracy of tba 
GPS updatast tba noisa lavels of tba position and velocity as well as 
tba filter knowledge of these components were increased first in ve- 
locity Vy a factor of 4 and run as case 25. Tben tbe position com- 
ponents were increased tqr a factor of 4« restoring tbe position error 
to its original value* Ibis was run as ease 25. All other parameters 
were as previously specified. 

Tbe results show a definite impact on performance resulting from tbe 
bigber noise levels in tbe GPS position and velocity channels * 

However, tbe increased sise of tbe error in performance did not corre- 
spond directly with tbe increased iseasurment noise* 

Neither attitude error nor its convergence characteristic were signifi- 
cantly affected by these changes in measurement noise or tbe resulting 
increased spacecraft navigational error. Tbe projected error in posi- 
tion on tbe earth's surface will be increased by tbe additional space- 
craft navigation position uncertainty in the manner previously 
discussed. 

Effects of Clouds 

Cloud conditions were simulated by picking a random number from a uni- 
form population* Ibis number yields a cloud density by interpolation 
from a density table. Should tbe cloud density exceed 40Z, tbe system 
attempting a GCP update will default to one of tbe star trackers not 
occulted by Earth, moon, or sun. If tbe cloud density is between lOZ 
and 40Z, tbe GCP update is vtade but with degraded accuracy from 3 
meters (normal) standard deviation on tbe surface to 30 iseters standard 
deviation on tbe surface. 


Case 30 examined tbe effects of clouds using the standard mechanisation 
with a linear cloud density table from OZ to 90Z* All other parameters 
were nominal, with GCP and GPS update intervals of 80 and 5 seconds, 
respectively* 

Tbe results are not greatly different from those previously seen except 
in attitude* It appears that the convergence rate is bigber in tbe 
region from 0 to 300 seconds. This seeias to result from tbe default to 
tbe star trackers occurring in tbe region of 400 seconds. Tbe attitude 
update with star trackers is not coupled to tbe position, thereby 
resulting in faster convergence* 

Near 1100 seconds on tbe attitude error curve, there is some increased 
attitude error. This resulted from cloud cover between 10 to 40Z, 
forcing tbe bigber GCP error referred to. Tbe results did recouverge 
to acceptable values quickly* 




CCP Itecovry Witte and Without Star Tracker Backup 


Case 31 aiaulated tba effects of long periods of tins without GCP | 

updates* Here the landaark tracker was turned off and the systen 

allotmd to default to the star trackers for attitude update* Case 32 

is siailar to 31 although the star trackers were not available for 

update* The period of shutdown for both cases was approxisiately 1000 | 

seconds after the systen converged (2000 seconds) from its startup < 

transient, 2000 to 3000 seconds* All other systen paraneters were held | 

at noBunal with a GPS updates at 5 seconds and GCP updates when active 

at 80 seconds* s| 

In both cases the attitude error tends to grow in the period of no 
landmark tracker updates. However the bound on attitude error in case 
31 is held down bP the star tracker. Each case converges rapidly once ^ 

the landmark tracker updates are again laade active* 


C* OVERALL CONCLUSIONS 


The system as sisnilated met the requirements set down earlier for con- 
ditions where the GCP update intervals are between 20 and 80 seconds 
with GPS updates taken at about 10-second intervals* Although the 
bounds appear marginal (12 to 14 meters, Table IV-5), data are taken 
close to the starting time by comparison with some systems and the 
mission duration* In some cases part of the initial transient was in- 
cluded in estimation of the peak or 3o value for a given parameter* 
Therefore most of the results are felt to be pessimistic with respect 
to the size of the error bounds* These results also do not reflect the 
estimation of such gyro parameters as gyro bias, misalignment, scale 
factor error and nonorthogonality* With estimation and compensation 
for these parameters, the results should be further improved* 


D* RECOMMENDATIONS 


These studies took advantage of the fact that the reference body 
coordinate system is defined within the landmark tracker or GCP 
sensor* Therefore the misalignment of that sensor is by definition 
equal to zero* However, this shifts the error contribution of mis- 
alignment of other sensors if such misalignments are not estimated and 
compensated* The effects of such uncompensated terms as gyro bias may 
be seen in cases where long GCP update intervals are used. Further 
I studies should be undertaken to implement the estimation and compensa- 

^ tion of constant gyro parameters and possibly the landisark tracker 

I parameters of bias and misalignment should that sensor corrdinate frame 

i not be taken as a body reference. 
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Hi« g«n«rac«d by tb« Landsat nultiapactral acannar (M88) auffar 

from a variety of gaomatric diatortioaa. Uaara of Landaat iaaga prod- 
ucta raquira that the temporal reglatrationa batvaan aeanaa ramaln laaa 
than half a pixel. Bagiatration of inagaa raquiraa the maaauraaant of 
diaplacamanta through the uaa of GCP datactioni calculation of gaomat- 
ric diatortion coafficianta, and raaampling of the ran data to provide 
iiaagary with unifotm pixel apaeing. If any of thaaa functiona are in 
error, the final image will have aoma ragiatration error. Thera are no 
procaduraa for determining the quality of ragiatration provided by the 
Goddard Space Flight Canter. Such a technique muat be developed to 
periodically amaaure the performance of the image correction ayatem ao 
modela that compute the diatortion coefficienta may be updated aa the 
MSS and satellite parameters vary. Under this contract a program was 
designed to perform the function of performance measureomnt. 

The primary element in the image registration quality measureaient tech- 
nique had been previously developed by Martin Marietta. Although the 
programs were developed to analyze the accuracy of various correlation 
algorithms, a secondary output is the offset between the two images at 
each ground control point (GCP) location. A detailed description of 
these programs has been included in Volume III. 

A great deal of insight into the image registration error sources may 
be gained by plotting a grid pattern that illustrates the offset be- 
tween two images. Although this technique will not detect systematic 
registration errors, it provides a tool to study the random or periodic 
errors. 

When examining the registration quality of a number of scenes, it is 
necessary to establish a reference image to which all others may be 
compared. To create the diatortion grid, it must be asaumed that the 
reference scene is a uniformly sampled orthogonal image (note that this 
assumption precludes the smasurement of tbe registration accuracy in 
the reference scene but will provide for a relative accuracy oieaaure- 
ment of the other scenes analysed). The reference image may then be 
separated into a grid pattern (Fig. V-1) where each of tbe lines are 
uniformly spaced and orthogonal. 



- 4 — GCP Locations 


Figure V-1 

Referenoe Grid Pattern 
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Hi« dittortion of cii« ocofMc may ba aaaaurod by roglatorlng tb« imagory 
at aacb of tiia grid croaalnga in tha rafaranca wlcb the imagary in eiia 
aeanaa baing atudiad. Tha ragiatration vaetort wtoicb ia than juat a 
daacripeion of tha ralativ^ of facta* providaa a maaauTamant of tha lo-’ 
cal diatortiona* It ia important to raaliaa that a raeognisabla GCP 
will not ba praaant at aach grid croaaing* Tharafora GCPa nay ba aa> 
laetad from anywhara within tha aubaraa. Tha eorralation analyaia pro* 
gram haa baan awdifiad to allow diaplay of a uniform *:rthogonal grid 
pattam ovar tha rafaranca acana and allow tha aalaction of a GCP in 
aach aubaraa in tha grid. Tha apaeing of tha grid pattam will nomi* 
nally ba 200 pixala ao 400 GCP ragiatrationa will ba idaally provided* 
A viaual control logic rapraaentation (VCLR) of thia procaaa ia ahown 
in Figura V**2* 
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Figure V~2 GCP Selection 
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Th« correlation analyaia program haa also bean modified to automatical- 
ly regiater the 400 CCPa from the reference acene on the correaponding 
image being atudied. For each CCP detection that doea not encounter a 
falae lock, a regiatration vector depicting the reletive offaeta of the 
t«M> imagea will be diaplayed with the baae of the vector originating 
from the center of the CCP location in the reference grid (Fig. V-3). 

A VCLR of thia proceas ia ahown in Figure V-4. 



CCP Regiatration 
Vector a 


Figure V~5 

Regiatration Veator Diaplay 


A VCLR of thia proceaa ia ahown in Figure V-4. 
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Figure V-d Diatortion Meaaurment Frogrm 
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It is interesting to note thst since OOSALA provides s meesure of the 
registration offset at a single point, distortion nasaureiDent nay be 
performed without any program modifications. However, the entire proc 
ess of computing the distort i<^ at each grid point and the distortion 
vectors must be done by hand* 
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VI 


CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 


One of the key bottlenecks associated with NASA's end-tO‘>end data 
management problem is the process of image correction. Resampling not 
only requires a tremendous amount of data processing but also requires 
the integration of volumes of data from every aspect of the mission, 

• including gyro data, tracking data, and imagery. 

With the development of the MLA, autonomous real-time image correction 
will become feasible. The distortions associated with the multi- 

* spectral scanner's nonlinear mirror sweep and nonuniform sweep period 
will be eliminated, leaving a well-defined distortion caused by earth 
curvature and look angle. The tremendous increases in computational 
capabilities both on the ground and aboard the spacecraft will also 
simplify the problem. 

A solution to the problem may be realised through either ground pro- 
cessing or spaceborne processing. However, even though ground-based 
processing may alleviate the immediate problem, in the long range adap- 
tive systems using onboard intelligence will be required. Therefore it 
may be advantage'^us to begin the transition to onboard automated 
systems now. 

Many benefits may be realized through onboard spacecraft control. By 
isolating the science sensor from the attitude limit cycle of the ve- 
hicle using a pointing mount, a curved focal plane MLA that provides 
uniform sampling of the image can be incorporated, thus reducing the 
need for along-scan resampling. This approach, although attractive for 
missions such as Landsnt where the revisit cycle is periodic, may not 
be attractive for missions requiring deterministic data acquisition 
through sensor pointing because a curved focal plane array boresight 
must always point at nadir. However, isolation of the sensor from the 
attitude li'nit cycle will simplify the image correction process even 
^ without the curved focal plane array. 

One of the keys to fixture remote sensing missions is the onboard deter- 
mination of sensor boresight position in earth-fixed coordinates. Both 
image correction and deterministic data acquistion systems will require 
this knowledge to be useful. Analysis has shown that inertial systems 
alone do not provide suffirient accuracy for several reasons: 

1) Ic is difficult to resolve misalignments between inertial sensor 
frames and science sensc<r frames; 

2) Misalignment between inertial and local earth frames are difficult 
to compute; 

3} Star trackers are not accurate enough. 
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Ground control point (GCP) sightings can be used to directly solve for 
these misalignments and hence reduce the error in boresight 
determination. 

Analysis has shown that a registration accuracy of 0.1 pixels can be 
achieved using an SSDA correlation algorithm. It is also possible to 
detect situations where the correlator locks onto a false target. 

These factors make the correlator an attractive navigation sensor for 
an onboard image correction system. 

A system comprising two NASA standard star trackers, a NASA standard 
gyro package, a GPS receiver, and the landmark tracker with a 3U-meter 
resolution was simulated. Analysis results indicate that with a GPS 
update frequency of between 5 and 10 seconds and a GCP sighting fre- 
quency of between 20 and 80 seconds, the position of the sensor bore- 
sight can be determined to within 15 meters. With the one exception of 
the MLA, the entire system can be implemented with existing tech- 
nology. The most restricting feature of the system is the tremendous 
storage requirements for the GCP data base. Assuming that GCPs may 
only ^ located on land (accounting for 20Z of the earth's surface) and 
assuming a GCP update frequency every 40 seconds or approximately every 
1 1/3 scenes, the total storage requirements may be found by: 

3 2 2 

Total earth surface area « 4 tt(6 x 10 )' km 

“ 4.5 X 10®, 

Total number of available scenes “ 4.5 x 10® km^ x 0.2 t 175^ km^ 

- 2.9 X 10®, 

Total number of GCPs - 2.9 x iO® x 3/4 

= 2 X 10^, 

Total memory required “ 2 x 10® GCPs x 32® pixels x 8 bits = 18 x 10® bits. 

GCP pixel 


Although this is a tremendous storage requirement, it is still within 
the limits of current technology. This volume may also be reduced by a 
factor of two by reducing the number of gray scales in the GCP refer- 
ence from 256 to 16. If histogram equalization is used to reduce the 
number of gray scales, the accuracy of i.e correlation may not be af- 
fected. This possibility should be investigated. 

It is estimated that processing requirements can be met with the use of 
two general-purpose flight computers and a hardwired correlator. The 
separation of functions would be directly analogous to the navigation 
and the registration processors discussed in Chapter II. 
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RECOMMENDATIONS 


If the landmark tracker is to be realised as an operational system, a 
«ell*planned development schedule is required* The critical tasks that 
must be performed include: 

1) Further analysis of the correlator is required to determine its 
sensitivity to seasonal variations. Several passes over the saom 
scene are required and the GCP analysis process must be repeated 
for each pass. It is also necessary to determine the types of 
landmarks that produce the best correlation. This infonaation will 
be invaluable in the selection of GCPs for an operational system; 

2) Research must be conducted to determine methods of reducing the 
storage requirements of the GCP data base. The most promising ap- 
proach seems to be a mixture of image enhancement and gray scale 
reduction; 

3) A study must be conducted to establish the configuration of the on- 
board processing network; 

4) The entire system should be breadboarded and ground-tested to re- 
solve peculiarities of the design; 

5) The breadboard system may be flomi on a low-cost aircraft flight 
test to more fully demonstrate feasibility; 

6) The breadboard system could also be incorporated into the IAS 
ground demonstration as an additional step toward a total system 
configuration. 

The final development stage of the system would be incorporation of the 
design into an advanced remote sensing mission such as OERS. 

The curved focal plane MLA is another concept that should be investi- 
gated. If feasible, the concept would dramatically reduce the proces- 
sing requirements associated with remote sensing missions. This inves- 
tigation can be separated into two tasks. First, the feasibility of 
the concept must be analyzed. Technical considerations include: 

1) Precision required in machining of the curved surface; 

2) Effect of variations in local earth curvature on registration 
accuracy ; 

3) Pointing requirements needed to ensure 1/2 pixel registration. 

The second task that amst be performed is a cost/benefit tradeoff be- 
tween deterministic data acquisition through sensor pointing versus 
image correction using the curved focal plane array. The two ap- 
proaches are mutually exclusive* 
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